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DETAILED DESCRIPTION OF TFCHNICAL PROGRESS 
Task 1. Salt Selection 
Introduction 

Storage of thermal energy at temperatures above 315 °C (600*F) may be 
desirable for a number of applications, such as — 

a. Conventional base-load power plants 

b. Solar Thermal Power Systems 

c. Cyclic- or batch-type Industrial pror :sses 

d. Vehicular propulsion using heat engines. 

Generally shaft power from a thermal engine and/or thermal energy itself 
can be produced by heat exchange with the thermal-storage media (Figure 1). 



A76030*^96 

Figure 1 . THERMAL STORAGE SYSTEM FOR SHAFT POWER APPLICATION 

Consideration of the temperature of the thermal input, thermal stability, 
compatablllty of the storage media with their means of containment and heat 
exchange, and of the desired end product (shaft power or thermal energy) 
generally sets the choice of storage temperature. For example, in the case 
in which the thermal input to the storage media, comes from either application 
a or b above, a useful temperature range of" the storage media is expected 
to be A50® to 535®C (850® to 1000®F). The production of shaft power is the 
highest quality use of such stored thermal energy; therefore, maximizing 
the storage capacity within a narrow temperature band in this range is thermo- 
dynamically desirable b<^cause thermal engine efficiency is a function of the 
temperature of the input energy (Figure 2). If the quality (temperature) of 
the thermal input to the heat engine is allowed to vary, costly englne- 
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generation control equipment would be required to maintain a constant power 
output from the plant and the overall plant efficiency would suffer. The 
temperature of the thermal input to the heat engine should remain constant 
to maximize its efficiency. 

Table 1 shows some choices for thermal storage media in the 450° to 535°C 
(850* to 1000*F) temperature range. The higher storage capacity of a phase- 
change material (PCM) for which the latent heat provides the major share of 
the total capacity is readily seen. Prom the table we observe that large 
storage volume iflll be required if sensible heat is stored in high-pressure 
steam. A typical refractory material such as will require 3 times the 
volume and 4 times the weight of a typical PCM. 



Figure 2. tfOBlt EFFICIENCY USING HEAT AVAILABLE AT TEMPERATURE T AND 

REJECTING HEAT AT 100*F 
(Adapted From Hausz and Meyer’ ° ) 

Compared with the other media, molten carbonate salt PCM have the 

V 

following potential advantages: 

a. The melting points are in the desired temperature range (450* to 535*C). 

This means that the relatively high latent heats of melting >135 Btu/lb 
(>350 J/g) can be utilized, leading to compact system sizes, and that 
the heat is available at a constant temperature (at the melting point), 
thus providing good quality input into a thermal engine. 
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b. ^ Their good thermal conductivities relative to other salts result in 

efficient heat transfer and thermal laanagement* 

c. Their low vapor pressures result in low-pressure design and potentially 
long usable lifetimes. 

Table 1. WEIGHT AND VOLUME OF MEDIUM REQUIRED FOR STORAGE OF 10* Btu 
THERMAL ENERGY [Storage Medium Operating Between 



450* and 535*C (850* 

and 1000*F)1 



Pressure 

Height , 
Ib/lO* Btu 

Volume, 
ft*/10® Btu 

Molten Salts 
(Carbonates) 

AtiDOspheric 

7500 

53 

Ste^ 

100 psia 

850 

7,000 

MgO 

Atmospheric 

22,000 

100 


Sigh- temperature, molten salt TES systems have received limited attention, 
and data dn^ their heat transfer characteristics, cycle life, corrosion, and 
^It stability are lacking. The program at the Institute of Gas Technology 
(IGT) is directed at obtaining some of this information for carbonate salt 
mixtures suitable for use at temperatures ranging from 450" to 535"C (850* 
to 1000*F). The purpose of the salt selection task (Task 1) was twofold: 

1. To identify promising salt systems by comparing thermophysical and cost 
characteristics of salt mixtures, *which have a melting point between 
450* and 535"C (850* and 1000*F). 

2. To select a model salt system for initial experimental work. The 
methodology and data generated for this model system can then be easily 
extended to other systems. 

Selection Criteria 

Numerous salts and mixtures of salts can be considered for thermal energy 
storage in the 450* to 535*C (850* to 1000*F) temperature range. If coat, 
availability, stability, nontoxicity, and corrosiveness are the prime salt 
selection criteria, then the available choices are limited to chlorides, 
bromides, hydroxides, fluorides, and carbonates. In further comparing the 
inorganic salts, numerous other factors must also be considered.^ These 
include the following characteristics. 

Thermal Properties 

The important thermal properties of storage materials include the heat 

of fusion, thermal conductivity, and heat capacity. These factors determine 

6 
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the suitability of a material for heat-of-fusion thermal energy storage on 
the basis of its storage capacity and heat-^transfer characteristics. The 
heat of fusion determines the storage capacity of a material . The thermal 
conductivity of its liquid and solid phases will determine the charging and 
discharging rates obtairable, and will affect the complexity (and cost) of the 
heat -exchanger design. Materials with high heats of fusion and thermal 
conductivities will be desirable as they represent high-capacity systems that 
can be charged and discharged rapidly without the need for complex and costly 
heat-exchange systems. High heat capacities, although of secondary Importance, 
are also desirable because they are indicative of storage capacity above the 
heat-of-fuslon and reserve capabilities of the system. 

Physical Properties 

Storagu^material physical properties of interest are density, volume 
change with temperature (volumetric expansion), volume change on fusion, 
viscosity of the liquid phase, and vapor pressure. Density and heat of fusion 
determine the energy storage density. Volumetric thermal expansion must be 
considered when designing the container to compensate for the storage-material 
volume change in heating from ambient to the operating temperr ' ire range. 

Host materials increase in volume upon fusion. If these volume changes are 
large, difficulties in containment may arise or voids may form in the storage 
material during solidification. This can have detrimental effects on the 
conducting properties, causing losses in efficiency, charging rates, and 
possibly capacity. Viscosity of the liquid phase ncir the melting temperature 
is important because of its effect on convective currents. These currents, 
brought about by temperature variations in the PCM, aid in heat transfer and 
in the homogenization of the storage material. The vapor pressure reflects 
the stability of the storage material and affects the design of the contain- 
ment vessel as well as the losses that might result from a leak in the vessel. 
Ideally, the heat-of-fusion storage material will have a high energy density, 
little volumetric t lermal expansion or volume change when fused, and low 
viscosity and vapor pressure in the working temperature range. 

Cost, Containment, Safety, and Stability 

The storage material cost, although only a small part of the overall 
system cost, must be considered because of its effect on the economic feasi- 
bility of a system. The toxicity and compatibility of the storage material 

7 

INSTITUTE OF CAS TECHNOLOGY 









3/78 


8981 


with the heat-transfer Dedlum must also be considered because possible system 
failures may release storage materials into surrounding areas or bring storage 
and heat-transfer materials into contact. Toxic fumes or violent chemical 
reactions resulting from such systoa failures can produce harmful effects on 
life and the environment. Containment of storage materials becomes a problem 
because of Che possible corrosiveness of molten storage materials aud tr>e 
volume changes accompanying temperature changes and fusion. If any of these 
cliaracteristics occur to a great extent, complex containment designs may be 
required. Storage materials that show little corrosive activity tdien used 
with low-cost containment materials are desirable to maximize the life of a 
storage system. Lifetimes of 20 to 30 years will be required if storage 
systems are to function economically in commercial systems. This lifetime 
requirement necessitates good stability of thermal and physical properties 
in a storage material. 

Eutectic Vs. tfoneutectic Mixtures 

Noneutectic compositions melt incongruenCly, i.e., the solid forms a 
liquid and another solid of a different composition. This effect may be 
advantageous in that the second solid phase formed in the liquid could act 
as a nucleating agent during solidification of the liquid, preventing super- 
cooling effects. However, because of different compositions, the solid and 
liquid will also have different densities; this may result in separation by 
settling of the solid, producing a wide composition variation throughout the 
storage medium. This problem could be alleviated by some type of agitation, 
but that would only lead to increased system costs. Eutectic compositions 
melt congruently, and consequently, this problem would not be encountered if 
one were used. 

Change in Melting Point as a Function of Oomposition 

A change in the melting temperature of a salt mixture as a function of 
its composition within the system is also of primary importance. Compositions 
within a storage-material system may vary as a result of corrosion reactions 
or selective vaporization losses. If the melting temperature increases 
sufficiently, fusion may not be obtainable with the heat source for tfhich the 
storage ayscem was designed to operate. In this case, storage capacity would 
be decreased significantly, and the system would actually be converted to a 
sensible-heat storage system. 
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Conpatlblllty With Volume- Change Additives and Conductivity Promoters 

Thermal conductivities and volume changes occurring during operation are 
important considerations in determining the practicality of using a salt as 
a PCM. Although these properties are characteristic of the salts, they may 
be modified in the storage system by the addition of inert materials to the 
salts. Because metals have relatively high thermal conductivities, the addition 
of metallic screens, vires, rods, or tubes to the salt could enhance its 
overall heat-transfer properties, although some of its capacity vould be lost. 
The volume change that occurs when the total storage medium is fused may be 
controlled by the addition of a volume-change control additive to the salt. 
Compatibility of the PCM with such additives should also be considered when 
selecting the heat-of-fusion storage materials. 

Candidate Salts 

Mixtures of 31 candidate salts are presented in Table 2 %/ith the experi- 
mental or estimated values of properties relevant to their applicability as 
heat-of-fuslon TES materials. These salts were selected for consideration 
because their melting points fall in or near the 450° to 53S°C (850° to 1000°P) 
range and they do not display any particular difficulties in handling, contain- 
ment, stability, or availability. Some salts that can be used for thermal 
energy storage at temperatures either well below or above the range of primary 
Interest have also been Included for comparative purposes. The materials 
under consideration include single salts and binary and ternary salt mixtures 
of both eutectic and noneutectic compositions. Quaternary systems were not 
Investigated due to the lack of pertinent thermal behavior data and suitable 
procedures for estimating the properties needed for evaluation. 

Cost ($/10* Btu) as a function of volume (ft*/10* Btu) of the salts 
presented In Table 2 Is shown graphically in Figure 3. The practicality of 
using very hlgh-capaclty, high-cost salts for large-scale applications Is 
questionable because of the Importance of salt cost In determining the degree 
to which a TES system can be economically competitive with alternative energy 
storage systems. The highest cost salts are composed of bromides, fluorides, 
hydroxides, or some chlorides (mixtures 1, 2, IS, 18, 21, and 25). The lowest 
cost salts are composed entirely of mixtures of alkali and alkaline earth 
chloride salts (mixtures 6 through 9, 19, 20, 22, and 23), %fhereas moderate- 
to-high cost salts are composed of carbonate salt mixtures (nos. 11 through 
14, 24, 28, 29, and 30), mixtures of carbonates, and some chloride and fluoride 
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salts (mixtures 5, 10, 26, and 27). Several of these lov-to-moderate cost 
salts emit highly toxic fumes when heated (mixtures 5, 6, 8, 19, 20, 22, 23, 

26, and 27), presenting potential safety hazards that would restrict their 
use and increase the cost of containment to prevent leakage. 

The heat capacity of the solid phase, reflective of sensible-heat storage 
capacity, and the thermal conductivity of the liquid phase, indicative of 
heat transfer rates, are plotted in Figure 4 for chloride and carbonate salt 
compositions. Carbonate salt mixtures are superior to chlorides in both of 
these properties by a factor of 2 in each case. 

Hixtures of lithium, potassium, and sodium carbonates or of sodium, 
calcium, and magnesium chlorides are the most practical salts fcr use as PCM. 
Hoiiever, the thermal properties of mixtures of carbonate salts, particularly 
their superior thermal conductivities compared with those of other salts and 
their smaller volume changes on fusion, suke them best suited as model storage 
materials in the 450* Co S35*C (850* to 1000*F) range. Mixture 28, a 35 weight 
percent Li2COj-65 weight percent ^ 2^3 ■ole percent Li2CX>^-50 mole percent 
mixture, %«s selected as a model system for experlsKntal work. This 
is a congruently melting mixture having a heat of fusion of 148 ^tu/lb 
(345 J/g), chat forms an intermediate compound, LiKCO^. Because experimental 
values of the thermophysical properties of the Li 2 C 0 ^-K 2 C 0 ^ mixture are well- 
characterized (compared with, for example, mixtures 12 and 29, which are less 
expensive) and because the change in melting point as a function of composition 
is very sauill, this mixture is Ideally suited as a model system. This mixture 
will, therefore, be used to define heat-transfer characteristice and potential 
problems and also to provide Che first-cut engineering data required for 
large-scale system design. 

Task 2. Heat-Transfer Analysis and System Design 

Introduction 

Heat-transfer analysis of TES systems using PCM is necessary to determine 
the interrelationships among the various factors affecting energy storage, 
to estimate the relative siagnitudes of various heat-transfer resistances and 
the overall heat-transfer flux, and to develop the mathematics and methodology 
required Co aid in the design of phase-change TES systems. After such an 
analysis is developed and is proven workable in the laboratory, accurate 
design and cost estimates can be made to determine the feasibility of storing 
thersul energy in PCM. 12 
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A heat-transfer analysis of a PCM storage system presents two major 
complexities. First, heat is stored In a PCM in both its latent and sensible 
forms, and the solid-liquid interface moves durin}^ the transfer of heat* 
Secondly, because of the regenerative nature of the storage system, a cyclic 
solution has to be considered. Due to these complexities, only a limited 
number of approximate solutions have been obtained in the past. We have used 
tnese solutions, with suitable modifications, to develop a heat-transfer 
model. In the following sections, we will discuss a) the hack;;round and a 
brief review ot past work, b) observed interri^lat itvish i ps ann)n^ variables 
of interest, c) experimental verification of seme approximat ions, and d) the 
design of an experimental unit using the approximate solutions and the heat- 
transfer model developed from them. 

Background and Past Wor k 

Basic Equations in P h ase- Ch ange Heat T ransfer 

A regenerative heat-storage and retrieval system is comprised of a finite 
mass of storage material and possibly two working fluids. One of these fluids 
supplies some of its heat content to storage, and the other fluid transfers 
the stared thermal energy from the lES system. In our ^^ase, the storage 
material undergoes a phase change, so that both the liquid and solid phases 
are present in the unit and the solid-liquid interface moves during heat 
transfer. Because the heat capacitier* and conductivities of the two phases 
are different and the boundarv between the two phases is not fixed in space, 
the mathematical analysis of hert transfer requires simultaneous consideration 
of both phases with an appropriate boundary condition Thus, heat flow In the 
solid phase follows Fourier*s law of heat conduction - 


* div grad ORIGINilL PACK U) 

s OP POOR QUAIJW 

In the absence of convective currents, heat flow in the liquid phase can be 
described by — 



div grad 


( 2 ) 


At the Bovlng intertace, the heat Is absorbed or released In the form of latent 
heat of fusion, so that - 


^ «rad t^ Rrad - .A 




O) 
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and 

t, - 1, - 1_ <*) 

The heat tranaferred to the working fluid ia giwan hy — 

<V» "f ^ <'f - V 

which la coupled by the relation — 

k grad t « U (t^ - t^) (6) 

at the wall of the heat-’tranafer coiMiuit. Solutions of the ayataa of 
Equationa 1, 2» and 6 exiat only for a few selected geometries, with additional 
si^plificationa introduced to make the solution process tractable. 

Additional complexities arc involwed becauae of the periodic nature of 
the storage^retrieval cycle. If the same fluid is used in the haatins and 
cooling half-cycles, the inlet temperature, t. can be repreoMted by a 

f f MM 

periodic function, such as — 

^£, in “ **hot ■*’ ‘cold ^‘bot “ ‘cold' 

where — 


‘hot 

‘cold 


nexlaum temperature of the fiot stream 
mlnlmnn temperature of the cold stream 
frequency. 


for a more general case, a Fourier series eumiiing the individual harmonica of 
Equation 7 can be used to construct the actual temperature wariatiM in t^ ^ 

The available analytical solutions can be claeaiftad in turn eategerieat 
1) tranaient (aoncyclic) solutions that asaome as li^tlally uniform tuaper- 
ature profile and 2) cyclic solutloaa that account for the reitamaratlve netura 
of the etoraga syatam. 


The transient solutions ate useful for experimantal studies because 
a) the heat-tnnefer coefficients can he determlnad or verified by recording 
the time-dependent temperature behavior of the wortUng fluids and b) heating 
the energy-storage materials past tlw> fdiaae chaaga into the li^id ^ase 
results in starting the solldificatiim half-cycles at approximately uniform 
tamperatures. He have found experimentally, and the hydro d yn amic st^iility 
theory also suggeets, that a temperature gradient ia the liquid phase of the 
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reference salt caK;ioc be aaintained in a layer that is more chan a fraccion 
of an inch Chick. This means Chat Che solidification halT-cycle can be 
studied analytically using cransieitL solutions for several geomecrlcal 
conf iguraC ions . 

Periodic solutions are i^>ortant for storage system design, but solutions 
involving phase changes are unknown in the literature. Useful Information can 
be obtained, however, using the available aolutions for regenerative heat 
transfer in senalble-bcat storage meterlals. 

Tranalent Solutions 

An approximate solution to two^phase beat transfer can be obtained if 
sensible heat ia assumed to be negligible compered with latent beat and If 
perfact contact is assumed between the coolant and the solid storage neterial. 
For a flat-plate g( 


try (Figure 5), approximating grad c by (t - t 1/a, we have 

ft % 


vdT ■ pAH^ da 




t ) At/s 

3 


and, upon Integrating — 


kT(c 

O.i e* * “ 




dAH, 


or In d iswnsionless form ~ 


0.5 !l 


ph 


(8> 


(9) 


( 10 ) 


whare — 


• Fourier number 


k ^ 

CO a* 
P 


AH, 


ph 


■ pbaee-cksngc mari>er * 
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In a hollow cylinder with R • r/r , where r is the r^tdius of the cooling 

o o 

tube, a siailar analysis leads to — 

■ro ■ *•* V 

■ 1 1 la cooling froa outside; R 2 1 In cooling fron inside. Bramlette *»t al .^ 

used these approrlaarr aolntions to coepare heat fluxes anong dlfierent 

I inif ainasnr configurations and inferred that an annular arrangeoent (working 

fluid inside) would provide the best heat fluxes. The assumptions Involved 

in Cbe heat balance, nanely c_ 0 and U • are not general. A check 

P 

agalnat a solution (to be discussed below) that does not resort to these 
asstn^fons showed that as calculated fron Equation 11, was off by a 

factor of 2 to 3 for values of between 9.5 and 10, using 
Ite do not recoanend using these approxinatlons in designing t storage unit. 


An analytical aolution for finite heat capacity, Cp, and AH with a 
flat-plate geometry was presented by Neumann,'* who assumed that a) the PCM 
was Initially at the nelting point., t|g, b) there was no temperature gradient 
in the liquid phase; and c) the surface at x • 0 was subsequently held at 
tn ^ A nr/re realistic model must include a finite heat-transfer coeffi- 
cient, h, at the boundary and conductance, k, across the wall; that is, 

U * Ho assumptions should be aade on the nature of the temperature 

profile in the solid phase; i.e. , a finite sensible heat of the solid phase 
should be included. Such an analysis was performed by Megerlin.^^ In his 
so! tion, the time-dependent trial functions appear as coefficients In a 
powi series expansion of the depth coordinate. 


Equation 12 gives the dimsnalonless time needed to solidify the molten 
material from r* • 1 to r* • R 
R 


r*ln r*(H^ln r* + 2) dr* 


r* D* + 2g^ In r* (R^^ln r* + 2) 


ORIGINAL PAG® Jfc 
CfP POOR QUALITY 


(Ng^ln r* + 1) (12) 


Solved for the solidified layer thickness R, which is measured in terms of the 
internal radius r^. Equation 12 yields the thickness of the solid-phase layer 
as a function of tine. Two Illustrative results obtained by numerical evalu- 
ation of the integral in Equation 12 a: / shown in Figure 6. The family of 
curves with ^laaa number llp|| « 0.57 describes the solid- layer thickness of 
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the reference naterial cooled with roos- temperature alr.*^ Phase-change nuad>er 
« 10 refers to the same material, except that the coolant runs about 
25^C colder than the melting temperature of the salt. 

These curves can be used to calculate various parameters, such as time 
required to solidify a given core of PCM, effect of heat-transfer coefficient, 
coolant temperature, and conductivity promoters. For example. If the coolant 
air is only approximately 25*C cooler than the melting point of the PCM 
(Np}j * 10), the time required for solidifying a given layer of salt would be 
approximately 10 times longer than that required if room-temperature air were 
used (Npii “ 0.57). On the other hand if time is kept fixed at Np^ * 10 for 
exaiq>le, the thickness of the layer r* should be reduced by half. Similarly, 
the effect of the relative resistance of heat transfer on the heat-transfer- 
fluid side and the storage-medlnm side can be determined by observing the 
effect of the Biot number * Uro/kg) on the thickness of the solidified 
layer. As can be observed from Figure 6 for Np|^ « 10 and a fixed thickness of 
a given PCM, an increase in from 0.5 to 5 results in a twofold decrease 

in the time required, T. This suggests that the heat-transfer resistance in 
the solid storage material is limiting for > 5. For -Z 0.5 (the value 
expected in our engineering-scale iinit), the heat-transfer resistance in the 
PCM still appears to be the major factor, although the contribution of heat 
transfer on the heat-transfer-medium side is also significant. If the con- 
ductivity, k, of the PCM is doubled, that is, if “ 0.25, the value of T 
required to solidify a layer of R >■ 6 decreases to 757, of the original value. 


Periodic Solutions 

As mentioned earlier, no periodic solutions are available for phase- 
change heat transfer. Therefore, we can consider a simple periodic solution 


tFor illustrative purposes, we used properties of the salt system LlKCO, 
selected for our experimental %n>rk. Relevant thermophysical properties of 
this system are*^^ 


Melting point : S05*C 

Heat of Fusion (AH^): 148 Btu/lb 


Heat capacity at melting point: (c„)g ■ 0.32 Btu/lb-*F 

(cj)i - 0.42 Btu/lb-*F 

Thermal conductivity at melting point: kg • 1.30 Btu/hr-*F-ft 

kg - 1.09 Btu/hr-*F-ft 


Viscosity at melting point plus 25**C: p « IS cP 

Density at 298*K; p ■ 141.5 Ib/ft* 

Density at melting point: «• 125.5 Ib/ft’ 
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for a flat-plate, sensible-heat system, and then introduce a concept of 
"apparent sensible heat" to account for latent heat and sensible heat 
simultaneouslv. 


For flat-plate geometry. Equation 1 can be written as — 


3t 


3^t 


3t “ “ 3x2 


(13) 


For a cyclic variation in the working-fluid temperature, the temperature in 
the PCM is given as — 


t(x, t) ■ Aq exp t”x(l + i) K) exp (wot) 

The thermal diffusivity a >= k/cpP can be modified by assuming the apparent 
Cp of a fictitious solid, which can be obtained as follows: 

Heat released per unit mass, Q, by a PCM between temperatures tg| and 
ta, having a heat capacity (cp)g, and latent heat AHf, is given as — 

Q “ (cp)s (^n, - t^) + AHf 

and for the fictitious solid with heat capacity (^p)app* ^ given as — 

Q 


(14) 


(15) 


(cp)app (tjn “ 


(16) 


From these two equations — 




(<=p>app * (‘^p^s f 1 + (cp)g (t„ - t«)- 


(17) 




or, in terms of N 


ph 


“ (C.)« (1 + 


(^p>app ■ ^^P's «ph^ 


(18) 


Note that the value of (cp)^pp is dependent on At * t^^ - t^ and is not purely 
a property of the PCM. Also, because during cooldown the entire mass of PCM 
does no., attain temperature t^, either an average temperature or the temperature 
of the working fluid can be used. When Np^ >> 1 — 

AH. 

Nph * 7— T- (19) 


^^p^app ^ (^p^s 




Temperature Response of a Storage Unit 
To estimate the temperature response of a storage unit undergoing cyclic 


operation, we used the periodic solution of Equation 14. The following values 
are calculated for our storage system: 

0.60 Btu/lb-®F (assumed; based on comparison with Neumann's 


^^p^app 


solution) 
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total period ** 24 hours for a daily cycle 
0.0153 ft^/hr 
3 ft“» 


8981 


app 


K 

X 


» 2.1 ft 


Phase velocity ■ Xv) ■ 0.0875 ft/hr “ 1 in. /hr 

This means that, in terms of order of magnitude only, the designer of a 
periodic energy storage system (daily cycle) utilizing LIKCO 3 as a storage 
material and operating between the icelting point (505*C) and ambient temperature 
(22*C) should think in terms of 2-foot-long toiperature waves penetrating the 
material with the speed of about 1 in. /hr. Note, however, that for any 
harmonic oscillation (for exaaqile, a sine-wave input), the agq>lltude of the 
incoming wave will be damped by exp (— 2 ir) • 0.0019 over a single wavelength. 

With such damping, choosing the thickness of the material as great as a whole 
wavelength would appear to represent underutilization and therefore waste. 

Thus, we would like to determine an optimum thickness or distance between the 
heat-exchange surfaces. 


In working with a slab of finite thickness 2a, it is convenient to define 
a dimensionless frequency M — 


M - 

M can also be expressed in terms of a Fourier 
thickness, a, and the total period — 



( 20 ) 

number based on the slab half- 


Therefore, M ■ 




OBIOmAL P4<?^ 
OF POOR QUALICy 


( 21 ) 

( 22 ) 


The temperature response of a slab to a cyclic heat input is usually 
expressed in terms of the parameter M. The relevant solutions are available 
in literature, but their numerical evaluation is laborious . Relevant 
solutions for the temperature response of a storage medium under periodic 
variation in a heat~ transfer fluid are presented elsewhere, ^ but an interest- 
ing observation can be made regarding the optimum thickness for maximum 
storage. Although very little heat can be stored in a thin slab ("a" or M 
small for a fixed period, t^), a semi-infinite body with its surface temper- 
ature varying harmonically can store, at most, 2 /tt (or 0.798) of the maximum. 

Groeber'’ observed that the heat stored during a half-cycle goes through a 
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mnximum of 0.911 at M * 1.2 and decr#»i.^*>s for V ^ 1.2 until the limit of 
0.798 is reached for practically all values of K > 3. This suftgests that, 
in the absence of other overriding considerations, the designer of a storage 
device with a flat-plate matrix might start at M « 1.2. However, tradeoffs 
between various competing parameters and economics should be considered. 

Experliaental Verification of Simplifying Assumptions 

A 5- in. -tall, 3-in.-diam laboratory model with a l/2-ln.-OD interior well 
was used to study solidification half-cycles with the model salt. For exact 
dimensions, constructions details, and thermocouple placement, see Figures 7 
and 8, This study confirmed two assumptions developed and used in the analy- 
tical treatment of the heat storage problem with a phase change: 1) Assume 
that there is no thermal gradient in the liquid phase, so that in a cyclic oper- 
ation the solidification half-cycle always begins with a uniform temperature 
distribution, permitting the use of transient (“single blow”) solutions, and 
2) use an “apparent” for system design. Both assumptions appear to be 
confirmed by the test results shewn in Figure 9. 

The No-Thermal-Gradient Assumption; First Part of the Cooling 
Half-Cycle (Cooling the Liquid) 

The model syster (L 1 KC 03 ) was heated to about 50**C above the melting point. 
Cooling begins with the liquid phase uniformly at S50°C. The cooling was 
provided by 22“C air passing through the 0.43-in. -ID internal well at 50 ft/s. 
Use of the Dittus-Boelter equation* for turbulent forced convection inside 


tubes in the form 


h - 0.20 (125) G*’* (0.43/12)'**^ 


where G is the mass velocity, G 


0.0742 (50) - 3.71 Ibm/ft^-s, and 


constant 125 is a temperature-dependent factor Involving air conductivity and 
viscosity (k/u ), gives a value for the heat-transfer coefficient of 13.9 
Btu/hr-f t^-*F. Several trial-and-error calculations of the exit air tempera- 
tures based on the relationship between the Increase of the enthalpy of air 
and the heat transferred from the hot charge along the interior well length z. 


P dz 


hp (t,^. - t) 


or, integrated. 


^exit *~w 
*-a ” *-w 


exp ( 
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Figure 7. SCHEMATIC OF LABOR.\TORY TES SETUP 
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lead to corrected values of h ■ 13 for the heat-transfer coefficient and 
65.6*C (ISO**?) = for the bulk air temperature. The knowledge of 

of the surface coefficient, h. allows us to calculate the te^>erature in the 
liquid phase using standard solutions for conduction heat transfer with boundary 
conditions of the third kind. A solution^ of these equations for hollow 
cylinders assuming a finite conductivity of the material (i.e. , no mixing) 
yields curves (Ai and Bi, Figure 9) that fail to reproduce the experimental 
results (curves A, B). 

The assumption of perfect internal nixing and nr temperature gradient 
leads to a cooling rate of 

_ i -8*F/min 5*C/min) (26) 

QT DU Cp tt 


which reproduces the experimental cooling rate of about S*C/min, and, on in- 
tegration to t * 150 *<* 872 exp (-0.0116 T), where T is in minutes, matches 
the experimental time- temperature curves A and B quite accurately. This 
suggests vigorous convective (buoyancy driven) mixing in the liquid phase in 
a cell of the size used here. 


In fact, the Rayleigh number “Assures the ratio of 

buoyant to thermal forces in a fluid. With the kinematic viscosity of the 
liquid estimated at V • 0.17 ft^/hr and the thermal diffuaivity ■ 0.0183 
ft^/hr, the Prandtl number v/a * 9 (comparable to that of water at room 
temperature). Taking At - S0*C (*90*F), the coefficient of volume expansion. 


8 - 


- P 


pAt 


1.4 X 10“'*R~‘ 


(27) 


Since the significant dimension L * 1/12, the Rayleigh number is then of the 
order of 10^. This is several orders of magnitude higher than the stability 
criterion (Ng^ ~ lO’). ^ 

To approach the condition of stability, the width of the liquid-filled gap 
would have to be reduced by a factor of 10, to about 1/10 inch. Subsequent 
experiments on solidifying the model storage system (LlKCOs) in a honeycomb 
matrix of about 1/4-in. mesh showed some temperature gradient between thermo- 
couples A and B but not to the extent predicted by the nonconvection solution 
curves Aj and Bj. 

26 OF POOR QUAUTF 


INSTITUTE 


TECHNOLOGY 



0 P 


GAS 



3/78 


8981 


Use of an "Apparent” Cpt Second Part of the Cooling Half-Cycle 
(The Phase Change) 

Figure ^ shows s brief (about 2.5-min) isothemal time lag at 500“C for 
thermocouple A and a 22-ain time lag for thermocouple B. Megerlin's solution 
(Equation 12, Figure 5) for Npjj = 0.57 and Ngj * 0.2 predicts longer solidi- 
fication times, 4.7 and 28 min, respectively. This solution, however, is 
written for an infinitely long cylinder with the heat sink at the interior 
well only. Vlhen heat losses from the outer lateral surface and from the top 
and bottom arc added to the central well loss and the cooling cates are re- 
calculated on that basis, the isothermal lag times become 2.7 min for thermo- 
couple A and 20 min for thermocouple B. 

Third Part of the Cooling Half-Cycle (C. oltng the Solid Phase) 

Thermocouple A begins to indicate solid phase at its location after 12.5 
minutes of cooling (Figure 9). The apparent Cp over a temperature range of 
505* to 445*C (941* to 833*F) is 0.32 148/108 =■ 1.69. Used in the analytical 

solution for the hollow cylinder, it results in the curve A-, whicii closely 
follows the experimental curve A over this range. The use of Cp * 0.32, i.e., 
of the sensible heat alone, results in Curve A';', which docs not agree with the 
experimental data. 

Thermocouple B is so close to the outer wall that by the time its position 
has changed from the liquid to the solid phase, the whole mass Is solid; there 
is no evolution of the larcut heat, and Cpg 0.32 used in the analytical 
solution adequately describes the behavior of thermocouple B (curve branch Bj, 
Figure 9). 

Design of the Engineering-Scale (S-kUhr^j^) Experimental Storage Unit 

The engineering-scale units were designed with the overall goals of the 
work in mind: 1) to verify the PCM cooling behavior predicted by the analytical 

solutions, particularly the movement of the solid-liquid interface as a function 
of discharge time, because of the commercial TES design parameters tltat can be 
determined from it (such as heat exchanger spacing), and 2) to accurately 
determine energy discharge rates, their improvement, and system efficiencies, 
all of which play a significant role in determining the feasibility of storing 
thermal energy in PCM. 


27 


INSTITUTE 


0 f 


CAS 


TECHNOLOGY 


3/78 


8961 


Based on the considerations discussed in the previous sections, particular- 
ly the results of Brasalette et al «.^ and the experiaentsl results froa the 
laboratory unit, we have selected eu euuula'. arraogeaent Cor the experiaental 
storage unit. This arrangeaent can be scaled up. A coaaercial slsed unit with 
a aultihexagonal cell arrangeaent can be visualized. In addition, this arrange- 
ment is the siaplest experimentally. The final unit design is shown in Figure 10. 
The internal radius of the central veil used for charging and discharging was 
chosen to be 1 inch so that the well could accoaaodatc an electric heating 
eleaent large enough to sustain the power levels and teaperatures required to 
charge the unit in a reasonable tiae period (10-14 hours at - 1 kUe) while aain- 
tainlng an annular area large enough to generate a fully developed turbulent 
flow of the coolant air without requiring high pressures. This radius is also 
convenient for coaparlng observed perforaance to predicted perforaance (r 0 • 1 ) > 
Electric heating was selected because it allows easy recording of the chargii^ 
energy and, therefore, accurate deteraination of the storage efficiency. 

The diaaeter of the heating rod used is 0.75 in., leaving an annular 
cross-sectional area in the well of 2.7 in.’ (0.0187 ft’). Fans were selected 
to move the coolant air at the rate of 4600 ft ’/hr, representing a velocity of 
(4600 ft’/hr) (0.0187 ft’)"’ (3600 s/hr)”‘ ■ 70 ft/s. The corresponding 
Reynolds nuaber is — 


„ vd (70) (1.25) 

**Re " V " (12) (16. 36 X 10"^) 


44,570 


(28) 


indicating fully turbulent flow. The heat-transfer coefficient under these 
conditions can be deteralned by calculating the Russel t nuaber: 


Hhu • • 0.023 Hge*** - 120.5 

f 

Therefore — 


(29) 


h • (KhuXI'j)/** 

• (120. 5) (0.015)/ (1.83/12) (30) 

• 11.9 Btu/hr-ft’-*F 


The corresponding Biot nuaber is ~ 

H. -ifo. 1U-.-9)- 

**B1 (12) (1.3) 

28 


0.75 


(31) 
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The analytical solution for a phase-change nuaber « 0.j7 (calculated 
froa the properties of the aodel salt); and an « 0.75 (Figure 6) was 
used to deterainc the container outer radius, alloving col^>lete solidification 
in a 6 to 8-hour period. A 6-hour discharge period represents a Fourier 
nuaber of — 


»Fo 



(0.0287 X6) 
( 1 / 12)2 


- 25 


(32) 


From Figure 6 for • 25, llgj • 0.75, and Rpj, ■ 0.57, a layer of salt 
approxlantely 4.2 in. thick (k~l) can be solidified in a 6-hour period. This 
represents an outer container radius of 4.2 1 > 5.2 in. Siailarly, an 8-hour 

discharge tins repreaents a Fourier nuaber of 33, a salt thickness of 4.7 in., 
auA a container outer radios of 5.7 in. (11.5-in. diaa). On the basis of 
these values anJ the availability of construction aaterials, a I2-in. container 
diaacter was selected (Sdiedule 40 12-in. -diaa pipe). It «ras desired to have 
1 ft’ of salt in the container at rooa teaperature, which would then fill the 
container to a hei^t of 15.9 in. Allowing a lOZ salt voluae inwrease with 
teaperature rise tc 941*F (discussed in the next section), the aolten salt 
would fill the container to a height of 17.5 in. (disregarding the expansion 
of the container). Therefore, a container height of 18 in. (at room tei^>er- 
ature) was selected. 


Ten themoeouples were located in the salt by inserting 1/8-in. stain- 
less steel tubii^t (closed at one end) through the container bottoa. Eight 
of the theraocouplea ware located in the sane radial plane, five of these 
themoeouples at half depth (-0 in. frea the container bottoa), one located 
at the heat-transfer surface, and one each at radial positions (R-1) of 
1.6, 2.4, 3.2, and 4.8 in. These themoeouples were used to deteralne the 
aovenent of the solid-liquid interface with tine for conpari&'>n with the 
Interface aovenent predicted by the analytical solutions. The three other 
themoeouples were used to study axial teaperature response an/, edge effects 
not considered in the analytical solutions, 1 >n> other theraocouplea located 
in a perpendicular radial plane at (R-1) positions of 2.4 in. (one on each 
side of the heat exchanger) and at v^^rying levels throughout the entire salt 
depth were uaed to study the axial teaperature response in detail. The 
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temperature Information generated by the thermocouples was used to approximat 
the shape of the interface as the salt solidified. The containers were also 
equipped to measure the height of the molten salt while in operation. 

The insulation was selected from the following consif rations. Under 
steady-state conditions the temperature distribution in a cylinder is given by ^ 

S 37> ■ » 

Integrating twice and using the condition at the insulation boundary 

-k ^ » h2(t2 " ta) (34) 

leads to 

ti - t 2 ^ h 2 V 2 £2 (35) 

t2 ~ ^ins 


where the subscripts denote — 

1 ^ outside radius, storage unit 

2 = outside radius, insulation 

a ^ ambient 

ins ■ insulation 

It Is ^esired» for safety reasons, and to keep heat los<^es low, to Veep the out- 
side insulation face at 54®C (130®F) or below; then t) = 505®C (941 ^F), 
t2 * 54 *C (130*F), ta « 22*C (72*F), h * 1.5 Btu/hr-ft^-"F (combined radiation 
and free convection), and “ 0.08 Btu/hr-ft-*F (mineral wool) leads to an 

t 2 of about 12.5 in., or an ii -sulatlon thickness (t 2 " ri) of 6.5 in. 

Verification of Theraophyslcal Properties of LiKCOa 

The purpose of this task was to determine or verify certain theriuophysical 
properties of LlKCOj critical to the performance of this salt as a PCM for 
thermal enorgy storage. Values for the melting point and heat of fusion have 
been reported in the literature and verified by DTA, DSC, and electrical 
conductivity measurements as a function of temperature. The melting point was 
found to be 505*C (reported value' ' of 504.5 t 1*(’)* with a heat of fusion of 
15i.5 Btu/lb (reported value‘s of 148 Btu/lb). The volume change on fusion was 
determined from measured thermal expansion coefficients between room temperature 

11 

INSTITUTB 0* CAS TECHNOLOGY 



3/78 


8981 


and 475*C (887®F;. A volumetric expansion of 4.5% was calculated from the 
expansion coefficients, representing a solid density of 2.13 g/cm^ at 505*C. 

The density of liquid LiKCOa at 505'’C was estimated from data of Speddlng^^ 
to be 2.01 g/cm^, and a volume change on fusion wae determined from these 
density differences to be 6 %. Therefore, a 10.5% volume change should be 
allowed for in containment design. 

Supercooling (common in low-temperature systems) was not observed at 
cooling rates above 5®C/min, and thermal cycling (in the laboratory unit) and 
corrosion products contained in the salt did not alter the fusion-solidification 
characteristics. It was also found that thermal cycling had a stabilizing 
effect on the LIKCO 3 . When the initial salt was examined by DTA, two endo- 
thermic peaks were observed that represent the 485*C eutectic temperature and 
the melting point of LiKCOa (Figure 11), indicating a composition slightly 
rich in Li 2 C 03 . After thermal cycling (12 cycles), DTA studies could only 
detect the fusion of LiKCOs at - 505*’C (Figure 12), indicating that the excess 
LI 2 CO 3 was selectively lost by vaporization, corrosive reactions with the 
container, or some other mechanism. 

investigation of Conductivity Enhancement and Volume-Change Control Additives 

Alkali metal carbonates have relatively higher thermal conductivities 
(■“1 Btu/hr-ft-'F) than the various inorganic salts considered. However, noting 
that the heat-transfer area required for a given quantity of heat depends 
strongly on thermal conductivity (particularly in the solid region of the salt), 
the heat-transfer area and the resulting cost of the storage system can be min- 
imized if the conductivity of the salt medium can be Improved. We therefore 
investigated the use of high-conductivity materials as conductivity promoters. 
Considerations for the selecclon of suitable conductivity promoters include — 

1. Thermal conductivity and heat capacity 

2. Compatibility with carbonates 

3. Availability 

4. Suitability of fabrication in desirable forms 

5. Cost. 

After the initial considerations, different materials in various shapes 
were experimentally tested in our laboratory TES (3-in.-diam container) unit. 
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The feaiiiblllty of using LlAlOa as a volune-control naterial to randomly 
distribute voids resulting from volume changes as the LilXOs solid^irs was 
assessed by examining the storage capacity sacrificed by the addition of 
LiAlOa* determining the rate at which the L5JU.02 particles settle from the 
carbonate^aluminate paste, and examining the influence of LiAlOa on the dis> 
cha..ge heat flux and system efficiency. The results of these studies are 
discussed below. 

Investigation of Conductivity-Promoter Additives 

As previously discussed under "Transient Solutions," the heat-transfer 
resistance in the PCM appears to be the major TES performance-limiting factor. 
This effect, which results because the low thermal conductivity in the solid 
phase damps the heat flux through the solid PCM formiNl on the salt-side heat- 
transfer surface, is the primary drawback in using a passive heat exchanger In 
such a system. If the thermal conductivity of the solid PCM can be cost- 
effectively increased, performance improvements will result. 

The most attractive method of conductivity enhancement appears to be by 
the addition of higher conductivity materials to the salt in the form of 
screens, wool, foams, honeycomb matrix, or particles. Several materials were 
considered for this application. Their relevant properties are summarized 
in Table 3. Some corrosion tests that were conducted with graphite showed 
that this material may not be stable in the molten carbonate environment. 
Dispersed particles do not provide significant heat-transfer enhancement, 
according to Siegel,^ ^ so they were not tested. 

To obtain relative Improvement in heat-transfer rates with the addition 
of conductivity promoter, we tested some available materials in our laboratory 
TES system. Figures 7 and 8 show the constructional details, and the locations 
of thermocouples in the system. A 3-ln.-dlam, S-in.-long stainless steel tube 
with the bottom end sealed, containing 800 g LIKCO 3 , was heated in an electric 
furnace with a heat input of 9 A at 115 V (heating rate of approximately 
5*C/mln). After the carbonate reached the equilibrium temperature, the cooling 
air (SO £t/s) was fed fiom the bottom of the l/?-in.-0D center tube. Testing 
was conducted with LIKCO 3 alone and for LIKCO 3 with a> 3X by volume stainless 
steel screen attached to the center cooling tube, b) 3X by volume aluminum 
honeycomb matrix immersed in the salt, and c) 3% by volume stainless steel 
wool immersed in the salt. The results of these tests, summarized in Table 4, 
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Table 4. C(^ARISON OF HEAT-TRANSFER RATES TO 
COOLING AIR TAOM LIKCO 3 WITH OR WITHOUT ADDITIVES 
(Conductivity Promoters and Volume Control Additl/es) 



Air Temp. *C 

Air Mass 

Operating 

O ^ 

r\ C 

q, ^Btu/ 
hr-ft* 


Description 



Flow Rate, 
Ib/min 

Temp Range, 

•c 

Btu 

Btu 

Qj/Q, X % 

LIKCO, 

24 

97 

0. 224 

545-450 

413 

281 

9281 

68 

UKCQ} with ^ voX % 
of SS Screen 

2t> 

106 

0. 222 

547-450 

420 

313 

10.340 

75 

UKCQ, wl«i 3 vol % 
of Aluminum Honeycomb 27 

116 

0. 222 

542-450 

412 

335 

11,060 

81 

LIKCCH with 3 vol % 
of SS Wool 

25 

119 

0.224 

548-450 

420 

365 

12. 0*^5 

87 

LIKCO, with 3 vol K 
LiAlOx 

24 

97 

0. 224 

550-450 

424 

279 

- 

66 

, t «40 

* .. 

T(t| dt. 

where 

tr integrated temperature of 

exit air 





t « time 


« period of time. 


« m^Cp(e) AT + AH*m^ 

♦ m c (m) At, where Q «theoreticel heat that can be extracted from pilot TES unit in the 
. ^ operating temperature range, 450* to 550*C« 

•mate of aalt, 1. 76 tb of LiKCO) 
c^( a) « apecific heat of aalt ( aolid or liquid) 

At « operating temperature range , *C 
ah .beat of fuaion. Ug Btu/lb for LiKCO, 
m^ • maaa of conductivity promoter 

specific heat of ronductivity promoter. 

^ Q| *' actual hrat transferred to cooling air from pilot Tt S unit in 40 minutes after cooling started. 

^ q « heat fiua ^uouhnw tube surface area «*.5| is,.* . 
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show that although stainless steel has a low thermal conductivity its 
conduct iivty-enhancing effect is good if used as a wool. We then obtained 
a relatively fine aluminum wool and tested it at concentrations of 3, 10» 
and 14 volume percent. 

The results of these tests are shown in Figure 13. The LIKCOb alone shows 
a uniform heat flux with cycling, whereas the units containing the aluminum 
wool show an initially higher heat flux, which Increases further with continued 
cycling. As a result of this behavior, a quantitative dependence of the output 
heat flux on the conductivity-promoter concentration could not be established. 
Further cycling of the units that contain 10% and 14% aluminum wool is expected 
to result in steady-state performance (no further change of the heat flux with 
cycling), as seen in the unit containing only 3% aluminum wool. When all of 
the containers are operating under this equilibrated condition, the appropriate 
relationship can be derived. The additional cycling, however, has not Teen 
performable because of mechanical problems. 

The primary observation to be made, therefore, is that at each concentration 
of aluminum wool, an increase in the output heat flux was immediately attained 
and it continued to Improve with cycling. 

Two possible expanatlons for the dependence of the output heat flux on 
cycling are — 


1. More intimate contact between the frozen salt and the aluminum wool re- 
sulting from slight corrosive reactions on the aluminum (possibly the 
formation of a thin LIAIO 2 surface layer), and 

2. A radial movement of the aluminum wool toward the heat exchange surface, 
resulting from the contraction of the salt as it solidifies. 

It is difficult at this time to ascertain which of these mechanisms may be 
controlling and wTiat contribution each makes toward the observed beltavlor. 
However, there is evidence for each, and the lesser depend *nce of heat flux on 
cycling observed in the salt containing the least amount of aluminum wool would 
support either. 


One of the primary difficulties in using a wool-type conductivity promoter 
is obtaining good contact with the heat exchange surface. Therefore, an alter- 
native means of improving the salt conductivity by ensuring good contact between 
the conductivity promoter and the heat-transfer surface was also tested. A 
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flnned-tube heat exchanger was tabrlcated and tested on the laboratory-scale 
unit. The eight fins, which were welded to the heat-transfer surface to ensure 
good contact, displaced ''2.5Z of the salt volume and extended midway into the 
salt. The results obtained from this unit are also shown in Figure 13. The 
maximum heat flux obtained is about the same as that with 3 volume percent 
aluminum wool. This results shows that a high-conductivity material distributed 
evenly through the salt in random orientation is as effective as a lower- 
conductivity material with a more Intimate contact. However, the fabrication 
costs are a great deal higher for a finned-tube heat exchanger. The aluminum 
wool is, therefore, as effective in enhancing conductivity, but is more cost 
effective, and appears to be the most promising conductivity-enhancement concept. 

Investigation of Volume-Change Control Additives 

The containment of a PCM used as a TES medium is complicated by the 
differences in thermal expansion between the containment material and the 
PCM and by the volume change accompanying fusion of the PCM. Differences in 
thermophysical behavior between storage and containment materials can lead to 
more complex containment design, increased containment cost, and lower volumetric 
storage capacity (based on container size). Each of these effects can decrease 
both the practicality and efficiency of TES that utilizes the heat of fusion of 
a PCM. The thermophysical properties of a PCM alone cannot be modified or 
controlled. However, they can be substantially reduced by creating a composite 
storage material containing both the PCM and an inert volume-control material. 
This latter material is characterized by its stability, thermal expansion, and 
particle size and morphology. The material is selected so that its thermal 
expansion is considerably lower than that of the PCM, thereby producing a 
composite thermal expansion lower than that of the PCM alone. 

Furthermore, the particle size distribution and morphology must be such 
that loosely packed particles form a tight, continuous capillary network. 

Vfhen the PCM is molten, all of the void space within the volume-control matrix 
is filled with the liquid. As the PCM solidifies and undergoes the accompanying 
volumetric contraction, the capillary network serves to localize the resulting 
voids on a microscopic scale, within the capillaries. The net eff> ct is to 
distribute the void space (resulting from the PCM volume change on solidifica- 
tion) throughout the composite storage medium. Therefore, on a macroscopic 
scale, the volume change accompanying a phase change and the void space 
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required to acconmodate it are barely detectable. As a result, any void form- 
ation or breakaway from the heat-transfer surfaces can be minimized. A closed 
containment vessel filled with a composite storage medium would not be subjected 
to the stresses induced by volumetric expansions and contractions accompanying 
fusion and solidification of the PCM. Therefore, associated complications in 
containment design and structure are eliminated. On the other hand, the 
incorporation of a volume-control material into a storage medium is disadvan- 
tageous in that the storage capacity lost by the displacement of the PCM is 
only partially recovered as sensible heat storage in the volume-control material. 

We have selected LiAlOa as a volume-control material because — 

• The PCM being studied (hlKCOs) is in the lithium carbonate-potassium 
carbonate system; the chesdcal inertness of LIAIO 2 to components in this 
system has already been proven; 

e The coefficient of thermal expansion of LiAlOj is approximately one-third 
that of LIKCO 3 ; and 

• The ability of this material to develop the required capillary network 
has been displayed in past work at IGT.^^ 

The effects on storage capacity (Btu/lb) and specific capacity (Btu/ft^) 
of lithium alumlnate additions to LiKCO) are shown in Figure 14. Capacity 
values for the composite :^torage medium were calculated for a lOO'C cycling 
(505* ± 50*0 . Under these conditions, the storage capacity of pure LlKCOs 
is 224 Btu/lb (148 Btu/lb as the heat of fusion and 76 £tu/lb as sensible 
heat). It is estimated that LIAIO 2 must be added in excess of 30Z by weight 
to effectively control the voluaw change of LlKCOs on fusion. The storage 
capacity of the composite is reduced to 170 Btu/lb with L 1 A 102 added in this 
amount. This composite represents a 242 decrease in storage capacity and a 
17.52 decrease in specific capacity from that of LlKCOs alone. 

More important is the effect of LiAlOs additions on capacity cost. The 
capacity cost of LlKCOs operating under 100*C cycling is $1790/milllon Btu. 

A composite storage medium containing 30 weight percent LiAlOs and 70 weight 
percen LlKCOs operating under the same 100*C cycling tiould have a capacity 
cost of $6150/mlllion Btu, a 2442 increase. 

Another anticipated difficulty is maintaining a uniform distribution of 
the volume-control material throughout the composite. If the volume-control 
material settles out of the composite, the capillary network needed to control 
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Che voluae change accoapanylng fusion of Che PCM is descroyed. As Che addicive 
■aceriaJ seccles, che voluaeCric changes of che coaposite increase wich cycling. 
This nay cause undesirable voids and a safecy hazard if the design for sCorage 
concainer expansion assuses chac che coaposice concains an even discribucion 
of che addicive aaCerial. The parcicle size and densicy of che addicive aacer- 
ial and che viscosicy and densicy of Che liquid in which che particles are sus- 
pended deteraine che rate at which the particles will setcle through che liquid. 
L 1 A 102 has a density of 2.6 g/cm^ (primarily 8 aud y crystalline phases); 
aolten LlKCOs has a density of 2.0 g/ca’ and a viscosity of 30 to IS cP near 
the aelting point. 

An analysis was made of this sediaentation effect. The settling rate of 
LIAIO 2 through aolten LlKOOa as a function of LIAIO 2 particle size was 
estiaated by Stokes* law and is shown in Figure 15. Sediaentation tests were 
conducted on a coaposite mixture containing 83 weight percent LiKCOs and 17 
weight percent L1A102 and having a particle size distribution between 0.01 
and 0.2u diameter. The powder tns packed into Type 306 stainless steel 
columns 1-1/4 in. in diameter and 12 in. deep. The columns were then heated 
to 530*C and air quenched after 3(X), 700, and 1250 hr at these teaperatures. 
Samples were then taken from the uppermost and lowest centimeter depths, and 
che L1A102 concentration was determined. The results of this test are shown 
in Figure 16. The settling rate calculated from Stokes* law for the average 
particle size of the L1A102 jsed in this test was 3 aa/10(X) hr. The settling 
rate determined from this test was Sam/lOOO hr, indicstlng that sediaentation 
would be a minor problem in long-term operation. 

Besides the problems of high cost and reducing the storage capacity, 
addition of L1A102 to LlKCOj can also hinder the heat-transfer characteristics. 

A composite storage medium containing 3 pe<'cent by volume LiAlOa in LIKCO) 
was tested in the pilot TES system to scudy this effect. The results are 
included in Table 4 for comparison wirH pure LiKCOj, and LiKCOa with conducti- 
vity promoters added. Taking pure LlKt'Os as a standard for comparison shows 
the heat recovery (Q^/Qg) reduced from 08 Z to 66 Z by the addition of L 1 A 102 
at this level, and the heat flux reduced from 9280 to 9215 Btu/hr-ft^. Al- 
though these changes are small, they Indicate the decrease in overall system 
efficiency resulting frcai the addition of L 1 A 102 to the HKCO 3 . 
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The decrease in storage capacity and system efficiency and the increase 
in capacity cost resulting from LiAlOa additions represent disadvantages of 
the use of LIAIO 2 as a volume-control additive. Because the volume change of 
LiKCOs in our system is < 10%, the benefits to be derived from volume-change 
control additive may be less important than the accompanying disadvantages. 
However, if a high-volume-change system such as chlorides is used, the volume- 
control additive may become Important. 

Construction Materials 

The heat exchanger conduits, the construction materials fox ''ontalnment 
of the salt, and the conductivity promoters must withstand the corrosive attack, 
of the salt. As we briefly described In our Task I Report on Salt Selection,^** 
carbonates are relatively less corrosive than chlorides, fluorides, sulfates, 
and phosphates. Molten carbonates are used in fuel cells, gasification of 
coal, ^ cleanup of sulfur-containing gases, ^ and other applications. A 
reasonable amount of corrosion data on molten carbonates if, therefore, avail- 
able. Although SF3clfic data lor our PCM, LiKCOs, is scarce, data on other 
alkali-metal carbonate mixtures have been reported; and we believe that the 
corrosion behavior of LIKCO3 would be similar to that of the other alkali- 
metal carbonates. 

Corrosion Data 

Janz and coworkers ^ ^ studied corrosion of Types 304 and 347 stainless 

steels in the Li 2 CC 3 -Na 2 C 03 -K 2 C 03 ternary eutectic under submerged conditions. 
From their experiments, the authors concluded that both Types 304 and 347 
are passivated under submerged conditions In an oxidizing environment. A 
corrosion of JO mils/yr was reported for Type 347 stainless steel in static 
tests at 730°C. (See Table 5.) 

Davis and Klnnihrugh^ studied corrosion by the Li 2 C 03 -Na 2 C 03 eutectic 
mixture under fully and partially submerged conditions, and reported a greater 
corrosion attack under partially submerged conditions. This observation Is 
Important for our TES system because a partially submerged condition Is likely 
to occur at the llquld-gas interface at the top of the liquid LlKCOs. 
Di.fferences in the activity of the oxidizing species in the meniscus area 
and in the bulk liquid are responsible for the higher degree of corrosion. 
Experimental verification of the expected greater corrosion will be discussed 
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in the following section. Davis and Kliinlbrugh also noted that ferro- 
alumlnum alloys, such as Kanthal A, Kanthal Al, and Hoskins 815 were more 
resistant than the 300-serles stainless steels.^ (See Table 5.) 

Atomics International^ has employed molten carbonates for sulfur removal 
and as reaction media for In situ combustion. They have carried out detailed 
studies of corrosion by carbonates under static and dynamic conditions; their 
melts consisted of the ternary eutectic by Itself as well as with additions 
of 20% sulfide or 20% sulfite melts. Some results of their 50-hour static 
tests and 1500-hour dynamic corrosion tests at 500**C are summarized In 
Table 5. The Atomics International report contains additional data for various 
alloys and ceramic materials. Their results suggest that — 

1. Up to 600*C, Type 347 stainless steel appears to be satisfactory for 
containing ternary carbonates in the presence or absence of sulfur 
compounds. A 1-year dynamic test at 500**C showed a corrosion of only 
0.2 mlls/yr fo** :his stainless steel. 

2. At higher temperatures, containment of these melts requires high- 
chromlum alloys, ceramics, cermets, or a frozen skull of the salt. 

Corrosion data on various grades of graphites has also been reported 
by Atomics International (Table 5). The corrosion rate obtained from 50- 
hour exposure tests varied from 20 to 150 mlls/yr. Graphite Is attractive 
as a conductivity promoter because of Its high conductivity, but Its 
corrosion behavior for different grades of interest needs to be Investigated. 

Aluminum Is also a possible candidate for a conductivity promoter be- 
cause It possesses a high conductivity and can also be considered as a 
construction material. Aluminum is resistant to carbonates at temperatures^^ 
up to 600*‘C; at the Institute of Gas Technology, we have operated molten 
carbonate fuel cells at 500** to 600**C using aluminum hardware.^ ^ Noting 
that the melting point of aluminum is bbO^C, the corrosion resistance of 
aluminum appears remarkable. This Is possible because of the formation of 
a protective alumina layer on the surface, which in turn reacts with carbon- 
ates to form a highly protective layer of lithium alumina te. ^ 

The TES unit should be designed for periodic (dally) cycles, so data 
on corrosion of materials under thermal-cycling conditions will be required. 
Atomics International also studied the effect of thermal cycles on corrosion 
rates. If the passivated films separate from the substrate during thermal 
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cycles, the corrosion rate is expected to be greater than that found in 

cons tan t-^temperature tests. For Haynes 25 and Hastelloy G and X, the corrosion 

rate was found to be 5 to 20 times greater under thermal cycling. 

Stress-corrosion cracking and the effect of chlorides, sulfides, sulfites, 
and sulfates under various condi t ions have also been studied/ Sulfur compounds 
enhanced corrosion, but chloride compounds did not alter the corrosion rate 
markedly. Therefore, it ii. desirable to avoid sulfur-bearing impurities in 
commercial LiKCOa powder. Higher corrosion rates can be expected if stresses 
exist in the components. 

Corrosion Mechanisms and Thermodynamics 

Examination of the outer film of the stainless steels showed^ that a 
dense layer of LiCr02 was formed, which is protective up to 700*^C. This con- 
clusion, based on X-ray diffraction and fluorescence studies, is in disagree- 
ment with that of Janz and Conte, who reported a LiFeOa layer as the corrosion 
product. It is possible that, because the X-ray diffraction patterns of LiFeOa 
and LiCr02 are similar,^ Janz and Conte could not distinguish between LiFeOz 
and LiCrO^- The results of Atomics International in using X-ray fluorescence 
are therefore more reliable. Note that, in either case, some loss of Li 2 C 03 
can be expected from the initial Li 2 C 03 -K 2 C 03 mixture. Because of the re- 
latively large quantities of salt mixture compared with the amount of corrosion 
products formed, the loss of Li 2 C 03 will be insignificant. 

Experimental Observations 

ORBmAL PAGE & 

The corrosion of two materials was measured: 0P pQOR QUALITY 

1. Type 304 stainless steel tube used for LiAlOa settling studies 

2. Aluminum honeycomb in the laboratory TES unit. 

Type 304 stainless steel tube was used for determining settling rates 
in a LiKCOs + LiAlOa mixture at 530*C (Figure 16). Two metallographic sections 
were obtained from the tube heated for 1250 hours: 1) in the meniscus area at 

the top of the tube, and 2) in the middle of the tube, where the interior 
wall of the tube Is completely surrounded by the molten salt. Figure 17 
show a metallograph of the material exposed in the meniscus area, which shows 
a corrosion of 2 mils in 1250 hours. This is equivalent to 14 mils/yr if a 
linear corrosion behavior is assumed. However, corrosion behavior is generally 
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Figure 17. METALLOGRAPH OF TYPE 304 SS IN THE 
MENISCUS REGION EXPOSED AT 530“C FOR 1250 HOURS (180 X) 
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parabolic [corrosion (time)^^]; thus a long-term corrosion rate much lower 
than 14 mils/yr can be expected. Compared with the observed corrosion in the 
meniscus area, the corrosion in the interior region was negligible. The 
observations confirm the differences in corrosion rates reported by Davis and 
Kinnibrugh^' for submerged and partially submerged conditions. 

Corrosion of aluminum honeycomb used as a conductivity promoter was 
observed in a terminated laboratory TES test that underwent 5 cycles between 
550® C and room temperature. Because the total exposure time was less than 
100 hours, the results are only tentative and cannot be extrapolated to predict 
long-term corrosion. However, an average corrosion of approximately 0.1 mil 
was measured. 

In summary, 300-series stainless steels (e.g.. Types 304 and VH) appear 
satisfactory under submerge! conditions. A somewhat greater amount of 
corrosion was observed on Type 304 stainless steel in the meniscus region, 
but additional tests are necessary to obtain long-term corrosion rates in this 
r "gion. Aluminum is an acceptable material, provided it is used in temperature 
ranges compatible with the desired mechanical properties. 
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Task 3. Construct, Operate, and Evaluate the Performance of the Englneerlng- 
Scale (8-kWhr^j^) System 

System Construction 

Based on the design considerations discussed in conjunction with Figure 10, 
3 engineering-scale units were constructed of 316 stainless steel. A cross 
section of one of the units is shown in Figure 18. The container is a 
Schedule 40 12-in. -diam pipe 18 in. long. The annular central tube, used for 
heating and cooling purposes, is a Schedule 40 2-ln.-diam pipe approximately 
24 in. long with welded flanges on both ends. The bottom of the container is 
a 1/4-ln. -thick circular plate. Heli-arc welding was used to secure and seal 
the bottom plate to both pipes. The thermocouple wells were initially install«!d 
by drilling 1/8-in. holes in the container bottom at the appropriate locations, 
inserting the 1/8-in. tubes to the proper height, and welding around their 
base. However, this method was later found unsatisfactory due to leaks 
developing around the welds, and the method was modified. Direct welding of 
the thermocouple wells was avoided by using butt-weld tube fittings (GYROLOK), 
which were found to perform satisfactorily. The container lid and outer 
sealing ring were machined from 3/8-in. plate stock. Sealing of the center 
tube was accomplish'‘d with a flange tightened by eight machine screws. The 
outer tube was sealed with a two-bolt compression-type flange. Aluminum gaskets 
were used on both flange arrangements. Standard 2- in. MSS slip-on flanges 
were used to connect the unit to the TES system. A schematic of the TkS 
system, comprised of the TES unit, an air blower, the heating element, the 
instrumentation, and the controls, is shown in Figure 19. 

A three-station test stand was constructed, with cooling air supplied 
through 2-in. -diam piping manifolded to each station from a single source. 

Air flow was adjusted and regulated by a limiting orifice control valve placed 
in each of the coolant lines. The units were Insulated with 7 in. of J-M 
Cerafelt mineral fiber. The insulation was applied in 1/2-in. -thick layers. 
Aluminum foil was placed between the layers to reduce radiation losses. 

System Operation and Data Collection 

The units were filled with approximately 130 lb of salt before mounting 
and insulation. rfiHna accomplished by a two-step procedure. First, 
the container was filled with powdered LiKC 03 , placed in an oven, heated to 
S80°C in an air/C02 atmosphere, held at temperature for 16 hours and then 
cooled. This step was repeated three times, filling approximately 80Z of the 
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Figure 18. ENGINEERING-SCALE UNIT MOUNTED IN TESTING STATION 
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container voluoe. The remaining salt was added by melting it in alumina 
crucibles, which were subsequently emptied into the container. The filled 
container was heated to 580°C and then cooled to room temperature. 

The units were operated under varying cycle periods (a cycle consisting 
of charge, holding time, and discharge). The units were charged with a 
3/4— in.— diam H-T flrerod fixed in the center of the central tube. The he-'itlng 
rate was controlled with a programmable controller > and the power input 
measured with a time-base recording wattmeter. The units were heated to 530*C 
and held at this temperature for varying times depending on the cycle period. 

The units were discharged with ambient air. Flow rates were determined 
using an orifice meter, inclined manometer, and an upstream static-pressure 
tap. Humidity of the air was measured by means of wet- and dry-bulb 
temperatures. The initial condition for a discharge run was for all thermo- 
couples at salt mid-level to reach 530*C. At the termination of the holding 
period, power to the heating element was cut, air-temperature thermocouples 
were inserted where required, md ai- flow was directed to the station under 
operation. Air temperatures were recorded continuously on a two-pen strip- 
chart recorder, and salt temperatures were recorded on a 24-point strip-chart 
recorder throughout the cycle. Air flow rates were checked and recorded 
periodically throughout the discharge run and adjusted when necessary to 
maintain a constant mass flow rate of -160 kg/hr (=^350 Ib/hr). 

The discharge run was ended when the slowest-cooling thermocouple at 
salt mid-level reached a temperature of 482°C. VHien this condition was 
reached, air flow was shut off, air-temperature thermocouples removed, and 
charging begun. 

The charging sub-cycle was controlled from a thermocouple located on the 
salt-side heat-transfer surface at salt mid-level. The charging program was 
designed to bring this location to the holding temperature, S30‘'C, In under 
an hour and therafter maintain this temperature throughout the charging period. 
A constant-temperature heat source, similar to that expected in commercial 
application, could therefore be simulated, allowing qualitative evaluation of 


Toaperature profiles through the salt obtained during shakedotm runs indicated 
that total solidification was achieved at a slightly higher temperature, 

493*C, at the slowest cooling location. 
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charging performance. The heater temperature was also recorded during charging 
in orde: to estimate the heat flux required to maintain this steady-state 
condition. Comparison with other temperature and flow requirements could then 
be made if alternate heat sources (such as high- temperature, high-pressure 
steam) were used for charging. 


Performance Evaluation 


Criteria for Performance Evaluation 


Several factors must be considered in evaluating the performance of a 
TES system, each of which has a significf*nt Influence on the engineering and 
economic feasibility of latent-heat energy storage. Although standards for 
performance evaluation have been established for lower- temperature systems^’ 
(primarily for solar energy applications), standards are not yet available 
for high-temperature systems such as that being investigated here. As a 
result, the evaluation must be based on the factors most influential on the 
performance of a system in a specific application (such as utility load-leveling). 
Therefore, we have selected the following parameters to evaluate the performance 
of the TES system being studied: 

• Average discharge heat flux over a discharge run. 

• Variation in the discharge heat flux with time during a discharge run. 



• Heat extraction efficiency — the percentage of the heat available to Che 
transfer fluid for a given AT and salt mass that is extracted by the 
transfer fluid. 

• Storage efficiency — the percentage of the charging energy recovered by the 
transfer fluid during a discharge run. 

• System stability — variations in system and materials behavior resulting 
from thermal cycling. 

• Cyclability, the capability of the system to be charged and discharged in 
the time period required as dictated by the anticipated duty cycles. 

Although other factors related to system performance have been considered 
and will be discussed briefly, those mentioned above will provide a strong 
base for TES feasibility assessment and relative evaluation of design modifi- 
cations (such as conductivity enhancement or alternate heat exchanger designs). 
Analysis and evaluation of system performance will be concentrated primarily 
on the discharge mode, because the discharge performance appears to be the 
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determinant %fhen the feasibility of commercial implementation is considered. 

This conclusion is drawn because the greatest barrier toward cost-effective 
system designs hinges on the ability to extract the stored energy in a time 
period which* for most TBS applications of primary interest, will be signifi- 
cantly shorter than the tine available for charging. Because the charging 
times are co^>aratively long and the initial formation of the higher-conductivity 
liquid phase occurs closest to the heat source during charging, the thermal 
conductivity of the solid salt (which is rate controlling during discharging) 
is a lesB significant factor during charging. 

Evaluation of LiKCOj Without Conductivity Enhancement (Engineering 
Scale) 

Two engineering-scale TES units were operated with LlKCO^ as the PCM 
without any conductivity or heat exchange enhancnmnt. The primary goals in 
operating these units were ~ 

a. To establish and verify baseline performance characteristics for relative 
evaluation of alternate PCM, conductivity-promotion, and heat exchanger 
concepts. 

b. To provide experimental verification of the mathematical heat-transfer 
models developed in Task 2. 

c. To provide first-cut engineering data required for the design of prototype 
TES hardware. 

The two units were operated for a combined total of SO cycles representing 
over 2500 hours of operating time. Both units displayed comparable and 
reproducible performance, and their cooling behavior was in close agreement 
with that observed in the laboratoi*y TES system and predicted hy che heat- 
transfer model. 

Heat fluxes were calculated from experimental data with the equation — 

q ■ mcn At/A (36) 

where ~ ORIGLNAL PAGE Ib 

q ■ heat flux, kJ/hr-m^ (Btu/hr-ft ) POOR QUALITX 

m • air mass flow rate, kg/hr (Ib/hr) 

• heat capacity of the air, kJ/kg-*C (Btu/lb-®F) 

At ■ air-temperature differential (t * - t. ), *C (*F) 

ouc in 

2 2 

A - heat-transfer area, a (ft ) 
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The calculated heat flux is therefore directly proportional to the air- 
temperature differential. The locations at which Che air temperature was 
measured in Unit 1 are shown in Figure 18 as A and A*. The temperature- 
measurement locations included B and B* in Unit 2, because the heat-transfer 
area used ico calculate the heat flux in both units was defined as the portion 
of the central cooling tube in direct contact with the salt, i.e., length BB'* 
Average heat fluxes were calculated to be of the order 6000 Bcu/hr-ft\ using 
air- temperature measurosents taken at A and A*. Similar calculations using 
air temperatures measured at B and B’ in Unit 2 resulted in average heat 
fluxes of 4000 Btu/hr-ft^ . 

The difference in heat fluxes determined from air- temperature measurements 
made at A and A* and B and B* could result from several factors: 

a. Heat transfer occurring in Che extended portions of the cooling tube. The 
heat available for transfer in these regions is that which is stored as 
sensible heat in the insulation, or that which can be drawn from the unit 
by conduction through the container walls, top, and bottom. The thermal 
conductivities of the insulation and steel, and the cross-sectional areas 
of the unit walls and cooling tube could allow enough heat to be trans- 
ferred to the extended portions of the cooling tube to account for about 
t%K>-thlrds of the observed heat flux difference. 

b. Different radiation effects on thermocouples. Temperature differentials 
measured between B and B' are not influenced by radiation, because thermo- 
couples at these locations are affected equally by radiant heating. 

However, a thermocouple located at A* %iould indicate a temperature higher 
than the actual air temperature, due to radiant heating, whereas a thermo- 
couple located at A would not be affected by radiation, because it is 
shielded by the coolant piping. Therefore, higher temperature differential 
would be measured between A and A*. 

c. Inlet air temperature at B is measured with a Type K thermocouple 
(Chromel-Alumel ) . Type T thermocouples (copper-constantan) , are used at 
A, A', and B* and are more accurate at lower temperatures. A Type K 
thermocouple had to be used at B because it remains in position during 
heating. A Type T thermocouple is not used because of the high temperatures 
reached (700*C) in the cooling tube during the heating period. Thermo- 
couples located at A’ and B' are not subjected to this heat because 

they are removed during the heating period. 

Although each of the last two factors may only make a small difference in 
the temperatures recorded at the (different locations, their combined effect 
could make a substantial difference in the calculated heat flux. This sensi- 
tivity of heat flux is a result of the TES unit Heslfe,.i> The heat-transfer 
area in the heat exchao»er portion of the cooling tube (B B’) is small (• 1 ft^). 
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%ihich results in a small air-tenperature differential at the air flow rates 
being used. Temperature differences of 1"C measured between A and B, and B' 
and A* combined can result in up to a 20X difference in the calculated heat 
flux during the latter stages of the cooling sub-cycle. Thi;. effect is less 
pronounced (3Z to 5Z) during the early stages of the coolxng sub-cycle when 
the At is large, but compounded over the entire cooling period it can result 
in substantial differences in the heat flux calculation. 

Performance evaluation of Onlt 1 is based on air-temperature measurements 
made only at locations A and A*. Efficiency and heat flu't therefore 
determined by heat tra.tsfer occurring in both the heat exchanger (of primary 
interest) and in the extended portions of the cooling tube. 

The calculated heat flux will be based on the heat exchanger area (BB'), 
because the temperature distribution along the cooling tube rectlons away from 
the heat exchanger is not known. As a result, the heat flux reported for 
Unit 1 will be an apparent value and not the actual flux occurring in the Heat 
exchanger. The actual heat flux could not be determined for Unit 1 because 
the appropriate air temperature * ^formation was not available. However, both 
the apparent and actual heat fluxes were determined in Unit 2, allowing 
verification of the performance observed in Unit 1 and providing accurate 
heat-transfer information for the TES system. 

Unit 1 Perf*/rmance — Heat Flux and Efficiency 

The first engineering TES unit contained 130 lb of LiKCOj and <as cycled 
20 times. After the twentieth cycle, the unit was cooled to room temperature, 
the insulation removed, and the unit inspected. The outside surface of the 
container was oxidized in the typical fashion consistent with heating stainless 
steel to about 538*C i. air, but no evidence of salt loss by creepage was 
visible. Slight vapor loss was apparent as indicated by discoloration of the 
insulation at the junction of the sealing flanges. The unit was only partially 
sealed with an aluminum gasket, and the low vapor losses observed show a 

23 24 

significant advantage of this system over systems using fluorides * as 
phase-change materials, where hermetic sealing is required to prevent moisture 
absorption and vapor losses. A leak was also detected in the container bottom. 
The exact location of the leak could not be determinv*:^, but is strongly 
suspected to have occurred in the weld around the thermocouple I'ell located 
near the cooling tube. Approximately 3.5 lb of salt was lost from the container 
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by leakage. The leaks appear to have a self-sealing effect, i.e., as the 
salt spreads through the insulation it cools and solidifies, preventing more 
salt from escaping. 

The container was opened to inspect the inner surfaces, particularly in 
the meniscus region where a higher corrosion is expected, and to obtain a salt 
sample for analysis. The solidified salt upper surface is shown in Figure 20. 







PV8041009 


Figure 20. SOLIDIFIED SALT SURFACE Of TES CONTAINER NO. 1 

AFTER 10 CYCLES 

A cooling history of the salt can be approximated by examining the relative 

sizes of the salt crystals on the surface. Near the cooling tube the crystals 

are fine and numerous, and show a definite orientation with respect to the 

cooling tube. This type of growth pattern is the result of a very rapid 

cooling. Midway between the cooling tube r^nd outer wall, the crystals are 

large and randomly oriented, resulting from a slower cooling. Ne^r the outer 

edge, the crystals are again somewhat finer, but show no preferred orientation, 

indicating a cooling rate higher tb?an in the middle, but much slower than that 

ne*^r the cooling tube. These obsei tlons are in correspondence with the 

temperature profile through the salt as a function of cooling tim* shown in 

Figure 21. The flat portion of this temperature profile through th^ still- 

molten salt also confirms the strong convective mixing effects observed in 

the pilot TES units. it 
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Figure 21. TEMPERATURE PROFILE IN THE SALT AT HALF DEPTH 
DURING THE FIRST CYCLE 
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When a salt sample was removed, voids were found under the salt surface 
midway between the cooling tube and the outer wall of the container* It was 
therefore assumed that a solid **skin** formed on the surface of the molten salt 
due to heat losses through the Insulation. The salt level dropped due to 
leakage exposing the variable-depth thermocouple wells. It was also observed 
that the wells were bent, giving further evidence that a solid **skin” forms 
on the upper salt surface during the earlier stages of cooling. If such a 
skin did begin to form, it was expected that the solids would settle through 
the liquid due to the density differences between the two phases. One of the 
primary reasons incongruently melting compositions were not selected for study 
was this anticipated phase separation due to sedimentation. However, if this 
phase separation does not occur when incongruently melting compositions are 
used as the PCM, the design implications would be strong. This would allow a 
wider range of compositions within a salt system to be acceptable for TES 
applications. Thermal energy would then be available over a temperature range 
(between the solidus and liquidus temperatures for the composition being used), 
which would allow "tailoring** of the PCM to fit the application and duty cycle 
of a TES system. 

The temperature profiles obtained on cooling during cycles 1, 10, and 20 
are shown in Figures 21, 22, and 23. When allowances are made for differences 
in the initial salt temperatures, three distinct trends as a function of cycle 
become apparent: 

1. The average salt-side heat exchange surface temperature increases for the 
same cycle period. 

2. The variation in the salt-side heat exchange surface temperature decreases 
for the same cycle period, 

3. The salt temperatures at the two locations farthest from the heat exchange 
surface are substantially lowered for the same cycle period. 

Each of these trends results from improved heat transfer both across the heat 
exchange surface and through the fused salt. The improved heat transfer across 
the heat exchange surface is brought about by the formation of an oxide 
corrosion layer on the salt side during cycling. The bonding of the salt to 
this layer upon solidification is stronger and more Intimate than the bond 
developed between the salt and the uncorroded steel surface. The improved 
heat transfer through the salt is most likely a result of the elimination of 
randomly dispersed voids during solidification resulting from controlled cooling. 
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Figure 22. TEMPERATURE PROFILE IN THE SALT AT HALF DEPTH 
DURING THE TENTH CYCLE 
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Figure 23. TEMPERATURE PROFILE IN THE SALT AT HALF DEPTH 
DURING THE TWENTIETH CYCLE 
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The results of these changes in system thermal behavior are that — 

a. More thermal energy is extracted during a given cycle period, 

b. Heat is available at a more constant temperature, and 

c. Heat is available at a higher temperature. 

These results are manifested as overall higher heat fluxes during operation 
and can be seen qualitatively by examining the coolant differential temperature 
(t^^^ — during different cycles for a constant coolant flow rate. A 5% 
to 20% Increase in the temperature differential is seen between cycles 1 and 
10, but there is only slight variation between cycles 10, 15, and 20. We 
concluded from these results that system equilibration occurred during the 
first 10 cycles, representing approximately 250 hours of contact between the 
heat exchange surface and the molten salt. 

The performance observed during the tenth cycle, which is representative 
of Unit 1, is summarized in Table 6. 

Table 6. REPRESENTATIVE THERMAL PERFORMANCE OF TES UNIT 1 

Cycle 10 


2,820 

19,240 

^.740 

26,800 

2,100 

1.000 

29,900 

24,260 

5,640 

31,100 

9.11 

53,425 

15.64 

81 

78 

45 

6,100 
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The heat lost to the ambient during discharge corresponds closely to the 
power drawn during the holding period of the charging sub-cycle (time during 
i^lch the molten salt Is held at a constant temperature, 530*C). This input 
power demand results from heat loss from the charged unit. Distributing the 
heat loss (5640 Btu) evenly over the discharge period (5 hr, 20 min) yields 
an electrical equivalent of 310 watts. The observed power drawn during the 
holding period varied between 325 and 350 watts. 

The storage efficiency here is determined from two values of the energy 
input to the system. The first, termed the charging energy, consists of all 
energy input above the ambient heat loss level as determined by the power 
drawn during the holding period. This is shown schematically in Figure 24 as 
the shaded portion under the power input curve. The second, termed the total 
energy input, consists of the charging energy and the heat loss during both 
the charging and holding periods. This is shown in Figure 24 as the total 
area under the power input curve (shaded area + crosshatched area) . The 
storage efficiency based on the charging energy is 78Z, whereas that based on 
the total energy input is only 45%. The efficiency based on the total energy 
input is of primary Interest for commercial units, but the value derived here 
on that basis is not a valid approximation applicable to larger units. This 
results from the container size and geometry, and the cycle characteristics, 
as discussed below. 

To evaluate the efficiency of the engineering-scale unit on the basis 

of total energy Input Is misleading because the ratio of the total surface 

area of the insulated unit to the volume of the unit Is very large. The 

effect of scale and container geometry on the efficiency can readily be seen 

from the heat loss Information generated from this unit. The observed heat 

loss to the ambient during the 5.33-hr discharge period totaled 5640 Btu. The 

2 

total surface area of the Insulated container is 25.5 ft . This leads to a 

, 2 

heat loss rate of 42 Btu/hr-ft . Assuming this loss rate to remain constant 
for the Insulation and thickness used, regardless of container size, leads to 
a total losik during a 24-hr operating period of — 

ENERGY LOSS - (42 Btu/hr-ft^ ) (24 hr) (A) (37) 

where A * total surface area of the Insulated container. 
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The capacity of a storage unit can be determined from the volume of the PCMf 
Its density, and Its heat of fusion by the equation — 


CAPACITY = p AH.V 
o f s 


where ~ 


* density of the PCM 

AH^ > heat of fusion of the PCM 

V > volume of salt 
s 

If It Is further assumed that some fraction, f , of the container volume, V , 

s ’ c 

Is occupied by the molten PCM, and the remaining volume by heat exchangers 
and/or void space. Equation 38 becomes — 

CAPACITY - V p AH,f (39) 

CO f s 

The heat loss can then be determined as a percentage of the storage capacity 
by taking the ratio of Equation 37 to Equation 39: 


% LOSS 


ENERGY LOSS 
CAPACITY 


X 100 


(42 Btu/hr-ft’’)(24 hr) 

(V^)(p^)(AH )f 
CO t s 


(A) X 100 


Substituting appropriate values for the experimental PCM (LIKCO 3 ; 

128.2 Ib/ft^ , AHj “ 148 Btu/lb), Equation 40 becomes — 

% LOSS = (|^) (~) (41) 

s c 

It can be readily seen from Equation 41 that for any given f » the % loss 

8 

during a 24-hr period is directly proportional to the ratio of the surface area 
to the container volume. The value of V for the engineering-scale TES units 
is 0.83 with Vc ■ 1.18 ft^ and A » 25.5 tt . The ratio of (A/V«) » 21.6, and 


the % loss 


(|^)(21.6) - 138%. 


The effects of scale on the heat loss are shown in Figure 25. For this 

diagram, the heat loss during a 24-hr operating period was determined from 

Equation 41 for storage units with a height- to-radlus ratio equal to that of 

the engineering-scale unit tested here, but with varying radii and with 

V ■ 0.75. Doubling the size of the container (1-ft radius, 3-ft height) 
s 

decreases the proportion of capacity lost to about 40%. In the capacity rangr 

of commercial Interst (1-100 MWhr^, ), this value is reduced to below 10%. 
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Figure 25. HEAT LOSS DURING A 2A-hour OPERATING PERIOD FOR LARGER TES UNITS 
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As a result of the strong dependence of efficiency on scale, it vould be 
erroneous to Judge the performance of a larger-scale unit based on the efficiency 
found for the engineering-scale unit, if the total energy input is considered. 
Therefore, efficiencies of other units will be based on the charging energy 
only. Although this method of evaluation does not consider the heat loss 
during charging and holding periods, it will still give a conservative 
approximation of the storage efficiency because the high loss rate is maintained 
during discharging, thereby limiting the energy recoverable by the heat-transfer 
medium. 

The performance curve obtained from Run 1-20 is shown in Figure 26. One 

of the primary teasons heat-of-fusion TES is advantageous is that the stored 

energy is theoretically available at a constant temperature — the PCH melting 

point. As a result, the temperature differential between the PCH and the 

heat-transfer fluid would remain relatively constant; and if the mass flow rate 

and input temperature of the working fluid did not change significantly, the 

heat flux would remain at or near a constant value throughout the discharge 

period. The resulting performance curve would, therefore, appear flat. As 

can be seen from Figure 26, the heat flux decreases from 12,000 Btu/hr-ft^ 

2 

at 5 min, to 7000 Btu/hr-ft at 60 min discharge time. The rapid decrear-' 

results from the thermal conduction barrier Imposed as the si.lt solidifies on 

the heat-transfer surface and Increases in thickness. The sensitivity of the 

heat flux on the solid-salt thickness can be examined from Figure 26 and the 

* 

cooling data obtained at locations 2, 3, and 7 that are summarized in Table 7 
from Run 1-20. 

Table 7. SUMHARY OF COOLING BEHAVIOR AT LOCATIONS 2, 3, and 7 DURING RUN 1-20 


Location 

Salt 

Thickness, in. 

Time to 
Solidify, min 

Q, Btu/hr-ft^ 
(From Figure 26) 

3 

1.6 

126 

5700 

7 

2.4 

200 

4900 

2 

3.2 

320 

4100 


These locations are used because they represent 3 different radial locat uvi'- 
Isolated from "edge effects" such as heat loss to the ambient. 
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The heat flux is shown in Figure 27 as a function of solid-salt thickness, 

and shows a linear relationship. The slope of the line shows a decrease in 

2 

the heat flux of 1000 Btu/hr-ft per inch solid salt. Although the zelatlcnship 
holds true afcjir two hours of cooling, it does not appear accurate during the 
earlier stages of discharging. Extrapolating to time • 0 would indicate a 
heat flux of 7300 Btu/hr-ft^. The observed value, which is > 12,000 Btu/hr-ft^, 
indicates an even higher sensitivity during the initial cooling period. These 
results indicate strongly the need for conductivity enhancement or the develop- 
ment of active heat exchangers. If the performance curve cannot he flattened 
by one of these methods, steady-state heat transfer can only be achieved by 
varying the working fluid flow rates or entering temperature, which would lead 
to mere complex and expensive system design. 

After establishing basellue performanc characteristics, the first unit 
was operated above the design temperature to evaluate its influence on 
performance. Operating from a higher starting temperature may be advantageous 
because of the further utilization of the sensible heat storage capacity of 
molten salt and the subsequent Increase in system capacity. Broadening the 
system operating temperature range without altering performance would also 
make the system more compatible with variations in the duty cycle or the 
amount of energy available for storage. Run 1-15 was initiated from a starting 
salt temperature of 575*C, apprc :lmately 45‘C above the design temperature. 

The storage capacity was increased by about 3400 Btu, or about 12.5%. The 
performance curve obtained is shown in Figure 28, with the performance curve 
from Run 1-20 operated from the design temperature for comparison. Operating 
from the higher temperature resulted in an increase in the discharge heat flux 
during the first 85 minutes of discharging, after which the performance 
appears nearly identical to that when the system is operated from the design 
temperature . 

During the first 85 minutes most of the heat removed is sensible; therefore 
the solid-salt thickness Increases only slightly due to latent heat removal. 

As a result higher heat fluxes are obtained because of the larger driving 
for'..e due to the higher temperature differential (between the sail and the 
heat exchange fluid) and a lower resistance due to a thinner solidified salt 
layer on the salt side of the heat exchanger tube. 
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This hypothesis seems to be supported by the duration of the enhanced 
heat flux observed in Run 1-15. After 90 minutes of discharging^ the molten 
salt has been cooled to a uniform temperature of 505®C, the melting point. 

The heat flux enhancement disappears after about 85 minutes of discharging. 

The concept that sensible heat "displaces*’ latent heat and subsequently 

inhibits interfacial movement during cooling is also supported by these results. 

Both runs were terminated when nearly identical thermal conditions through 

the salt were reached. As can be seen in Figure 28, an additional 75 minutes 

was required in Run 1-15. Since the same amount of energy was extracted iti 

the form of latent heat and sensible heat removed from the solid, the difference 

in discharge times results from the different quantities of sensible heat 

stored in the liquid phases. If this "excess" sensible heat delays the removal 

of latent heat until the end of the cycle, it would be evidenced by a longer 

discharge time, which is calculable. The heat flux after 325 minutes of 

discharging (the time required to discharge the unit, starting from the design 

1 

temperature) is about 41C0 Btu/hr-ft , which corresponds to energy removal 
at the rate of — 

(4100 Btu/hr-ft^) (0.72 ft^) = 2950 Btu/hr (42) 

The additional 3400 Btu of excess sensible heat in the form of latent heat 
would require — 

(3400 Btu)/ (2950 Btu/hr) - 1.15 hr (43) 

or 69 min to be removed. This value agrees well with the 75 min observed. 

Unit i Stabil Itv and Cyclabilitv ^ 

1 1 JL OP POOR OUAIOT 

The first unit was in operation for a total of 1100 hours; and although 
this is a short period with respect to desired system lifetimes (> 44,000 hr), 
the performance stability observed Is encouraging. The effect on the average 
discharge heat flux of passing through the 20 cycles completed with this unit 
is shown in Figure 29. As observed in the laboratory unit (discussed under , 

Task 2), the heat flux appears to increase during the Initial cycling, but 
remains relatively constant after 7 cycles have been completed. This effect 
most likely results from enhanced contact between the salt and the salt-side 
heat-transfer surface brought about by slight oxidation of the heat-transfer s 
surface during oper.'':ion. After the completion of 7 cycles, the time required 
to discharge the unit also stabilized to an average of 5 hours and 40 minutes. 
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The discharge time varied by no more than 20 minutes during the remaining 
12 runs, when the system was operated from the design temperature. 

Chemical analysis, DTA, and DSC tests on the salt after the twentieth 
cycle showed no signs of degradation of the salt properties. The salt composi- 
tion, melting point, and heat of fusion remained constant during the operating 
period. It is interesting to note that the salt showed discoloration due to 
corrosion products after the salt filling and shakedown runs, but this was not 
evidenced in the salt removed from the unit after the twentieth cycle. This 
observation is important in that it indicates passivation of the container 
inner surfaces in a relatively short period of time. 

The system was readily cyclable in a 24-hour period. Full charging, 
accomplished in 8 hours, required an average charging flux of about 
5450 Btu/hr-ft^ (based on the charging power). A charging flux of about 
3800 Btu/hr-ft^ was required for a 12-hour charge. 

Increasing the holding time did not appear to affect the performance, 
indicating a good adaptability of the system to variations in duty cycles. 

Unit 2 Performance " Heat Flux and Efficiency 

Unit 2 is essentially a duplicate of the first unit, but incorporates 
improvements in temperature measurement and insulating techniques. The primary 
purpose for operation was to verify the system performance characteristics 
observed in the first unit and to study the heat exchanger performance in 
greater detail. The unit was operated in the same manner as Unit 1, but was 
discharged using higher coolant flow fates. A total of 30 cycles were completed, 
representing 1400 hours of operation, periods were varied by increasing 

holding and charging times, and the unit was discharged from both Initially 
fully charged and undercharged states. Stable performance comparable to that 
of Unit 1 was observed under all of the operating conditions. The system 
yielded average discharge fluxes of 6300 Btu/hr-ft with system efficiencies 
of 75% to 85% with air temperatures measured at A and A’ (Figure 18). However, 
when air temperatures were measured at B and B', the heat flux was found to be 
about 4000 Btu/hr-ft^ , which more accurately describes the heat exchanger 
performance. Strong convective mixing, which was also observed during dis- 
charging, enhanced the heat-transfer rates during the initial discharging 
periods. Detailed discussion of these results follows. 
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To verify the performance observed in the first er.^ineering-scale unit, 
several discharge runs were made on the second unit, with inlet and outlet 
coolant temperatures recorded at locations away from the system (A and A' in 
Figure 18) as was done with the first unit. Run 2-20 is representative of the 
observed performance and is summarized in Table 8. The designation indicates 
that Unit 2 was operated for 20 cycles. 

Table 8. "APPARENT" THERMAL PERFORMANCE OF TES UNIT 2 

Run 2-20 

Heat Flow, Btu 


Sensible Heat Removed From the Liquid 2,457 

Heat of Fusion Removed 18,944 

Sensible Heat Removed From the Solid 4,220 

Total Heat Extracted From the Salt 25,621 

Sensible Heat Removed From the Container 2,100 

Sensible Heat Removed From the Insulation 1,100 

Total Heat Lost by the System 28,821 

Heat Extracted by the Coolant Air 24 , 550 

Heat Lost to the Ambient 4,271 

Charging Energy, Btu 30,260 

kWhr 8.86 

Heat Extraction Efficiency, % 85 

Storage Efficiency, % 81 


The performance curves of Runs 1-10 and 2-20, based on air-temperature 
measurements made at A and A' in both units, are shown in Figure 30. The 
average heat fluxes were determined from graphical integration to be 
6300 Btu/hr-ft^ for Run 2-20 and 6080 Btu/hr-ft^ for Run 1-10. The difference 
between the heat fluxes for these runs can be accounted for by the difference 
in flow rates. Unit 2 was operated with an average air flow rate of 355 Ib/hr, 
approximately 7% higher than the flow rates used in Unit 1. Theory predicts 
a 4.5% increase in heat flux with a 7% Increase in the flow rate. The observed 
difference, 3.6%, agrees well with the predicted value and indicates that 
performance is affected by the flow rates being used for testing. The higher 
observed efficiencies result primarily from better insulation of the second 
unit, which decreases the heat loss during discharging. However, we concluded 
from the close agreement in heat flux between Units 1 and 2 that the observed 
performance is reproducible. 
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As previously mentioned, the performance discussed above is an ’’apparent** 
performance for the system. Although the values derived here are valid for 
comparison between the first two units, they do not represent the true per- 
formance, A better evaluation Is obtained by measuring the coolant temperatures 
the inlet and outlet of the heat exchanger (B and B* in Figure 18). 

The height of the molten salt was determined by an electrical conductivity 

method developed at IGT to verify the heat-transfer area used in the heat flux 

determination and to determine the precise locations for the outlet thermocouples. 

The measured salt height was found to agree within 1% of the value predicted 

26 

from density-temperature relationships and the thermal expansion properties 
of the containment materials, and the outlet thermocouples were located 
accordingly. 

Several runs from Unit 2 that were analyzed on this basis showed average 
heat fluxes varying over a narrow range between runs. The performance curves 
for two consecutive runs, which were made on a fully charged unit (2-15) and 
a partially charged unit (2-16), are shown in Figure 31. Run 2-15 results, 
representative of Unit 2 performance, are shown in Table 9. 

Table 9. REPRESENTATIVE THERMAL PERFORMANCE OF TES UNIT 2 

Run 2-15 

Heat Flow, Btu 


Sensible Heat Removed From the Liquid 2,312 

Heat of Fusion Removed 18,944 

Sensible Heat Removed From the Solid 3,727 

Total Heat Extracted From the Salt 24,983 

Sensible Heat Removed From the Container 2,000 

Sensible Heat Removed From the Insulation 1,050 

Total Heat Lost by the Syscem 28,033 

Heat Extracted by the Coolant Air 15,191 

Charging Energy, Btu 30,600 

kWlir 8.96 

Heat Extraction Efficiency (Heat Exchanger Only), % 54 

Storage Efficiency (Heat Exchanger Only),% 50 


The overall average heat flux for this run was determined by graphical 
Integration from Figure 31 to be 4024 Btu/hr-ft^'. 

The efficiency evaluation is complicated by the heat transfer that occurs 
In the extended portions of the cooling tubes. When an energy balance is made 
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on a run using air-temperature information collected at the heat exchanger 
inlet and cutlet (BB'), 40% to 50% of the stored heat is shovm to be not 
recovered by the heat- transfer fluid, and therefore lost to the ambient. 
Similar calculations made with air temperatures measured at the system in.’et 
and outlet (AA') show losses to the ambient to be half as large. These values 
correspond more closely to the heat losses predicted by input power drawn 
during holding periods. Although only 50% to 60% of the energy is recovered 
in the heat exchanger, it must be noted that about 70% of the stored energy 
not recovered in the heat exchanger is composed of heat transferred to the 
coolant outside of the heat exchanger, as evidenced in the first unit and 
shown in Run 2-20, Since no approximation can be made as to where this heat 
would have gone if the system were designed differently, it will not be 
considered in evaluating efficiency. As a result, the efficiencies will 
appear low, but should be used for comparative purposes only. 

An energy balance cannot be made on Run 2-16 because the exact charged 
capacity cannot be determined. A capacity value can be approximated, however, 
by comparing the total energy input during charging for Runs 2-lo and 2-16. 

The energy absorbed by the system during the changing sub-cycle of Run 2-15 
totaled 15.4 kWhr. This is the energy required to supply the sensible heats, 
heat of fusion, and heat losses to the ambient during the charging period. 

The energy absorbed in Run 2-16 during the charging period totaled 13.33 kWhr, 
or 86.5% of the energy drawn in Run 2-15. Assuming the charged capacity of 
Unit 2 in Run 2-16 to be 86.5% of that in Run 2-15, the following analysis 
can be made: 

Run 

2-16 


Heat Available for Extraction [ (0.865) (28,033) ] , Btu 24,248 

Heat Extracted by the Coolant Air, Btu ^ 13,326 

Heat Extraction Efficiency (Heat Exchanger Only), % 55 

Graphical integration of the performance curve for Run 2-16 shown in Figure 31 
yields an overall average heat flux of 4732 Btu/hr-ft“ . A similar measurement 
made for Run 2-15 over the same discharge period as for Run 2-16 (0-235 minutes) 
yields an average heat flux of 4382 Btu/hr-ft . We concluded from the close 
agreement in average heat fluxes and heat extraction efficiencies found between 
Runs 2-15 (fully charged) and 2-16 (partially charged) that operation of the 
system below rated capacity does not alter the system heac fluxes. 
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This is an important observation when the periodic nature of TES system 
operation is considered. Radial temperature profiles through the salt obtained 
during a 14-hour charging period are shown in Figure 32. The salt closest to 
the heat exchanger begins to thermally respond to charging within a half hour 
after the charging is initiated. However, the bulk of the salt continues to 
cool throughout the first 3 hours of charging. If in this case only 11 of the 
required 14 hours was available for charging, the system would only be charged 
to about 60% capacity. Although the discharging time would be decreased, 
stable performance would still be obtained. This factor would become even 
more important if demand dictated increased utilization of the sensible heat 
stored in the solidified salt because of the lowered salt temperature at the 
completion of the discharge. The subsequent charging time would, therefore, 
be Increased. 

One of the considerations made in selecting LlKCOa as an experimental PCM 
was the low liquid-phase viscosity at and above the melting point. This low 
viscosity enhances convective mixing in the molten salt, which improves the 
overall heat-transfer characteristics of the system. A strong mixing effect 
was indicated by the inability of a thermal gradient to be sustained through 
the molten salt as evidenced by the flat radial temperature profiles observed 
during discharge. The convective effects were also Indicated from axial 
temperature profiles ootained during discharging, shown in Figure 33. During 
the first half hour of discharging, the mean temperature in the upper two- 
thirds of the molten salt has decreased while the lower third has increased, 
after which a steadily decreasing flat profile is maintained to the melting 
point. The shape of the axial profiles also verifies the formation of a 
solid skin on the upper salt surface (discussed under "Unit 1 Performance — 

Heat Flux and Efficiency") and Indicates simultaneous solidification from 
the bottom and top. We anticipated the solid phase to form around the heat 
exchanger in a conical geometry based at the container bottom. However, the 
combined radial and axial profiles indicate an "hour glass" geometry that 
results from heat loss tlirougl' the top of the unit. 

Unit ? Stability and Cyclabillty 

Stable performance was observed in Unit 2 throughout the 1400 hours of 
operating time. The average heat fluxes observed during discharging as 
determined by air-temperature information collected from around the heat 
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Figure 32. SALT RESPONSE DURING CHAP.CING 
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Figure 33. VERTICAL TEMPERATCKE PROFILE AS A FUNCTION OF TIME AT A MIDRAD lAL PLANE 
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exchanger varied over a narrow range between 39o i and 4100 Btu/hr-ft^ when 
discharging was accomplished from and to design conditions. The discharge 
time required for each run Is shown In Figure 34. time required for 

discharge of the 15 runs, which started from design conditions, remained 
relatively constant, averaging 5 hours and 15 minutes, about 25 minutes faster 
than for the first unit, due to the higher flow rates. The charging performance 
was similar to that observed In the first unit; however, charging and holding 
times were varied over a wider range. No difficulties In completing a cycle 
were found if a minimum of 8 hours was allowed for charging. Charging In 
less than 8 hours required salt temperatures in excess of 593*’C in some locations 
during the charging period, thereby Increasing the possible corrosive reactions 
and vapor loss due to local overheating. 

All materials In the system also performed well ;hroughout the operating 
period. Corrosion was light on the interior of the unit, and no leaks were 
found. Slight vapor losses were indicated at the sealing flange junctions 
as In Unit 1, but no weight loss cou.'d be detected. The salt showed no signs 
of discoloration after 30 cycles, and chemical analysis showed no change In 
the salt composition within limits of detection (rO.5*). OTA and DSC testing 
of salt samples taken from six different locations In the unit Indicated well- 
defined fusion characteristics with melting points of 505” to 506”C and heats 
of fusion between 152 and 155 Btu/lb. No secondary compounds could be detected 
from thermal effects. 


Evaluation of LIKCO.. With a Finned Tube 


Unit 3 Operation — Flnned-Tube Heat Exchanger Design 

A third engineering-scale unit was constructed and operated as were the 
first two units, but It contained a flnned-tube heat exchanger to evaluate 
the effects of heat-transfer enhancement. The fins were constructed of 304 
stainless steel flat stock, 15 inches long, 3 inches wide, and 1/8 Inch thick. 
Eight grooves spaced every 45” were milled Into th«t cooling tube. The fins 
were placed in the grooves, and each was hell-aic welded around its entire 
base to ensure good heat transfer. They were spaced 3/4 inch from the bottom 
of the unit and 2 Inches from the top to avoid increased losses to the ambient. 
The relative radial spacing of the fins can be seen in Figure 35. They Increase 
the salt-side surface area of the cooling tube from about 100 to 705 square 
Inches . 
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Figure 35. FINNED-TUBE HEAT EXCHANGER LOCATED IN THE 
ENGINEERING-SCALE TES UNIT 

Unit 3 Performance — Heat Flux and Efficiency 

The third TES unit was cycled 11 times, representing 650 hours of operation. 
Performance evaluation made on this unit was based on air-temperature information 
collected at the heat exchanger inlet and outlet only (B-B* in Figure 18). 

The performance observed in Run 3-10 is summarized in Table. 10. 

The overall average heat flux for this run was determined by graphical 
integration to be 5250 Btu/hr-ft^ during discharging. 

The Influence on performance of the finned heat exchanger is summarized 
in Table 11, where average performance statistics for Units 2 and 3 are 
compared. The resulting performance Improvements are best explained by 
differences in the heat exchanger thermal state Induced by thtj addition of 
fins. The average salt-side heat exchanger surface temperature and surface 
temperature differential from top to bottom for Units 2 and 3 are shown as 
a function of discharge time in Figure 36. The use of fins (Unit 3) results 
in Improved conduction through the solid salt and produces a higher average 
surface temperature distributed more evenly along the entire length at the 
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Table 10. REPRESENTATIVE THERMAL PERFORMANCE OF TES UNIT 3 



Run 3-10 

Heat Flow, Btu 


Sensible Heat Removed From the Liquid 2,204 

Heat of Fusion Removed 18,944 

Sensible Heat Removed From the Solid 4,326 

Total Heat Extracted From the Salt 25,474 

Sensible Heat Removed From the Container 2,000 

Sensible Heat Removed From the Insulation 1,050 

Total Heat Lost by tne System 28,524 


8 

I 

I 

I 


Heat Extracted by the Coolant Air 17,170 

Charging Energy, Btu 28,938 

kWhr 8.47 

Heat Extraction Ef ffclency (Heat Exchanger Only) , % 60 

Storage Efficiency (Heat Exchanger Only), X 59 



Table 11. PERFORMANCE CHANGES RESULTING FROM ENGINEERING-SCALE 
OPERATION WITH HEAT-TRANSFER ENHANCEMENT 




Unit 2 

Unit 3 

Change 
With Fins 

Discharge 
Heat Flux 


4000 

Btu/ht-f t * 

5300 

Btu/hr-ft ^ 

+32% 

Heat Extraction 
"fficlency 


56% 

61% 

+9% 

Storage 

Efficiency 


53% 

60% 

+13% 

Discharge 

Time 


315 

min 

230 

min 

-27% 


heat exchanger. Because the heat-transfer rate Is directly dependent on the 
temperature difference between the heat-transfer surface and the coolant, a 
more efficient utilization of the heat-transfer area results; and subsequently 
a higher heat flux Is attained. The Improved thermal conduction resulting 
from the addition of the fins Isalso evidenced by a flattening of the radial 
temperature distribution (Figure 37) when compared to units without fins 
(Figure 38). The resulting Increase In the rate at which the solid-salt layer 
Increases In thickness Is shown In Figure 39. The Increased storage efficiency 
Is primarily a result of the shorter discharge time, decreasing heat losses. 
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Figure 37. SALT TBMPERATURE PROFILE AT HALF DEPTH DTOING 
DISCHARGING WITH A FINNED-TUBE HEAT EXCHANGER 
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The 32% improvement In the heat flux corresponds closely to the 25% Im- 
provement observed in the laboratory TES system (Figure 13 and discussion) when 
a finned- tube heat exchanger was used. As a conqparlson of the performance 
curves for Units 2 and 3» the averages of the last two runs made on each unit 
are presented in Figure 40. The goal of conduction enhancement is to both 
Increase and stabilize the heat flux during discharging. The performance in 
Unit 3 indicates the heat flux Improvement, but shows a more rapid change in 
the fltix with time after the bulk of the molten salt has been cooled to the 
melting point. This effect results from the diminishing Influence of the 
dimensions and spacing of the fins as the solid salt Increases in thickness, and 
is intensified by the shorter discharge time. 

The Influence of the fins is much less pronounced during charging. The 
radial temperature distribution during charging is shown for Unit 3 in 
Figure 41. Comparison with temperature profiles observed without any conduction 
enhancement (Figure 32) makes it evident that the fins Improve the charging 
characteristics only during the first 2 to 3 hours, or until enough molten 
salt has been formed to provide its own conduction enhancement. The net 
result is a 7% decrease in the time required. for charging. 

Unit 3 Stability and Cyclability 

Because of the short operating time, a quantitative assessment of tho 
unit stability was not practical. However, a qualitative indication of the 
system stability was found by the uniform discharge time per cycle shown in 
Figure 42. As previously discussed, both the charging and discharging times 
were reduced by the use of the finned-tube heat exchanger; and cycling in 
24 or 48-hour periods was achieved, although it required longer holding periods. 
This unit was not disassembled after completing operation; therefore materials 
performance was not established. 
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Figure 40. EFFECT OF FINS ON HEAT FLUX 
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Figure 41. SALT RESPONSE DURING CHARGING WITH FINNED-TUBE HEAT EXCHANGER 
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Task 4» Compare Experlnental Performance of Engineering-Scale System With 
Its Predicted Performance 

Introduction 

One of the primary goals of this program was to develop a heat>tir«-t>o; o 
model capable of predicting the thermal behavior of a TES system in < erms of 
Its method of operation and the properties of the storage material. Because 
phase change Is employed as the storage mechanism, a moving-boundary solution 
governed by heat-transfer characteristics was needed. Solutions of this type 
derived by various Investigators were examined for their thoroughness and 
applicability to the problem at hand as determined by the system under Invest- 
igation ^sre. (See section covering Task 2.) One such solution of phase- 
change heat flow proposed by Megerlin^^ was found to meet the criteria. In 
his solution, Megerlln derived an equation giving the dimensionless velocity 
of a solid-liquid Interface with respect to a heat-transfer surface in terms 
of the system design with suitable constraints. We converted this solution to 
a more usable form to give the dimensionless time (Np^) needed to solidify a 
molten material to a dimensionless thickness (R-1), and then evaluated It 
numerically for the conditions under which the experimental TES system was to 
be operated. The experimental system was dei^lgned to generate the Information 
required to verify Megerlin’s solution in our form. The results of this effort 
are discussed here. 


Solidification of LlKCOa Alone; Movement of the Solid-Liquid Interface 
During System Discharge 


An annular arrangement was selected for the unit design for ease of 
comparison with the theoretical model. Thermocouples were located at half 
depth In the salt at various radial locations (0.25, 1.6, 2.4, 3.2, and 4.8 
Inches from the heat exchanger), and their response with time was recorded 
during discharging. A typical discharge trace is shown with the unit design 
In Figure 43. The time required for the freezing front to reach the location 
of a thermocoup'.e was determined from the traces as the point where that 
thermocouple began to indicate a temperature below the melting point and sub- 
sequently yielded a descending cooling curve characteristic of a solid being 
cooled. The times required for the solid front to reach thermocouples isolated 
from edge effects (such as heat loss) are summarized In Table 12 for several 
standard runs made on Unit 2. The averages of these times are summarized In 
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Table 13 with the radial locations of the thermocouples and the calculated 
Fourier numbers. The theoretical solidification profile determined by numerica] 
evaluation of the Megerlin solution for our system (Npjj = 0.57, Ng^ = 0.75) is 
shown in Figure 44 with the experimental profile. The observed behavior is in 
good agreement with the predicted behavior during the early stages of discharge 
(Npo but is only fair during the remainder of the discharge period, which 

shows the solidification proceeding slower than predicted. 

The sources of deviation between the theoretical and experimental behaviors 
become apparent when the differences between the physical model and the 
idealized description given by the Megerlin solution are examined, as in 
Table 14. The primary source of erior appears to be in the excess enthalpy or 
superheat stored as sensible heat in the molten salt and container at the 
operating teiuperature (~20®C above t,„) . The superheat may be considered in the 
model two ways as follows. 

In the first treatment, the superheat is considered to affect the system 
behavior as if its source were the heat of fusion of the PCM. The sensible 
heat contained in the salt over the 20®C superheat amounts to ~17 Btu/lb. 

Adding this heat to the heat of fusion of LiKCOa results in an "apparent" 
heat of fusion of 165 Btu/lb (148 Btu/lb + 17 Btu/lb). Recalculating the phase- 
change number on this basis leads to an apparent Npj^ of 0.63, shifting the 
theoretical curve (Np|^ = 0.57) in Figure 44 to the right and yielding closer 
agreement with the experimental curve. If the heat content of the container 
above tj^ is included, the apparent Npj^ - 0.7, which brings the theoretical and 
experimental curves even closer together. 

In the second treatment, the concept of sensible heat displacing the heat 
of fusion is considered, as discussed under "Unit 1 Performance — Heat Flux 
and Efficiency." It has been calculated, from the knowledge of the average 
heat fluxes and the heat contents of the steel can and the insulation, that 
the loss of the specific superheat plus the loss of the can-and-insulation 
heat content in time, decoupled from the solidification process, amounts to 
about 34 minu. a, or 0,56 hour. If it is assumed that little or no solidifi- 
cation occurs until the superheat is removed, the solidification time scale is 
shifted by 0.56 hour, representing a correction on NpQ of 2.31 (Np^J = Np^ — 2.31). 
Applying this correction to the calculated Fourier numbers shown in Table 13 
(where applicable) almost brings the experimental and predicted curves in 
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Figure 44. COMPARISON OF SOLIDIFICATION PROFILES 
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Table lA. COMPARISON OF MEGERLIN SOLUTION 


Condition 
The Arrangement 

Internal Radius 

External Radius 

Axial Dimension 

Initial Temperature 

Excess Enthalpy 
Above tji 

Mlsc Losses 

Heat Content of 
Containers and Insul 


Mathematical Description 

Cylindrical, circular, 
euclsymmetrlc 

>^o 

Infinite 

Infinite 

At melting point tjg 


0 

0 

0 
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Physical Model 

Cylindrical} 
circular, axlsymmetrlc 

r^ » 1 In. 

Finite, r - 6 In. 

Finite, L “ 18 In. 

Above melting pt 

(0.42) (40) » 17 Btu/lb 
Via Insul, supports, etc. 

12 Btu/"F for the can 
22 Btu/^'F Insulation 
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Figure 44 together. Other variations, which. If corrected for, would tend to 
bring the curves even closer together are — 

a. The finite axial dimension of the physical model as opposed to the In- 
finitely long analytical solution suggests that, at the radius nearest 
to the entry of the coolant, the freezing front would form faster. (The 
coolant enters at 49“C, Is about 62*C in the middle at the thermocouple 
level, and leaves, on the average, at 77"C.) 

b. The convective currents in the liquid phase (left out of account completely 
In the analytical model) would tend to slow down the solidification process 

With these considerations in mind, the agreement between the observed and 
predicted behavior appears good throughout the entire solidification period. 

Use of the i!eat-Transfer Model in Large-Scale System Design 

Introduction 

The phase-change heat flow model developed by Megerlln^^ and evaluated at 
IGT during the course of this project can be used In large-scale system design. 
In the ICT form, Megerlln's solution determines the extent of thickness around 
a cylindrical Internal well (tube) carrying a cooling fluid to which a PCM will 
solidify as a function of the time it has been cooled from Its melting point, 
measured in terms of the radius o. the Internal well. The model accounts for 
the flow and thermal properties of the working fluid and the' thermophysical 
properties of the FCN (heat of fusion, heat capacity, melting point, thermal 
conductivity, and density). It yields the amount of a PCM that Is solidified 
with respect to the system design and mode of operation and, therefore, in- 
directly gives a measure of the energy transferred from the storage material 
to the working fluid under the same constraints. 

If the conceptual system Is a large molten salt bath with a number of pipes 
running parallel to each other in a well-formed matrix carrying the working 
fluid through the salt bath, the model will determine the pipe spacing, the 
number of pipes required, and their length (given the pipe radius) for a par- 
ticular storage requirement. Consequently, the model relates the energy dis- 
charged from a system to the physical design of the system, allowing a balance 
*;o be made between the energy required as determined by the working fluid 
conditions and the energy discharged that results from the system design. 
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TES System for Steam Upgrading 

To show how the heat-transfer model may be used in system design, we will 
apply it to a storage system in which high- temperature, high-pressure steam 
is th' working fluid. In this system, LIKCO3 will be used as the PCM for the 
purpose of superheating the steam at 1000 psi from 550*F (288“C) to 900®F (482®C) 
over a 6-hour period during which the steam flow will be 5.0 X 10® Ibm/hr. The 
enthalpy change in the steam is 249.2 Btu/lbm (from the Steam Tables), so the 
total system capacity must be — 

(249.2 Btu/lbm) (5 X 10®lbm/hr)(6 hr) » 7.48 X 10® Btu (44) 

which represents — 

(7.48 X 10® Btu)/ (148 Btu/lb) = 5.05 X 10® lb (45) 

of salt (heat of fusion only). The volume of salt that must be solidified 
(using the room-temperature density for a first approximation) is — 

(5.05 X 10® lb)/(141 Ib/ft®) * 3.58 X 10® ft® (46) 

This sets the "fixed volume" requirement that the heat-transfer model will be 
applied to solving. The phase-change number determined for a mean steam 
temperature of 725"F (385®C) is — 

« = ^ = ... - BtuAbJ 2 (47) 

Ph Cp(t^-t) (0.32 Btu/lb-"F) (941“F-725"F) 


The Megerlln plots for Np^ > 2 and a range of Blot numbers are shown in 
Figure 45. 

In this example, we have chosen to use pipes with a 2.5-ln. diameter, and 
need to determine the number of pipes required, their spacing, and their length. 
The Fourier number corresponding to a 1.25-in. (0.104-ft) radius and a 6-hr 
discharge time is — 

M , 24 , (Q-0287. ftVhrmtoJl . 3^3 g 

NFo ^ (0.104 ft)^ ^ ^ ^ 

Given the pipe diameter, 2r{,, the Blot number need only be specified now 

to use the Megerlln plots in Figure 45 to determine r,, which will yield the 
thickness of thesolldlfled salt (r^ - ro) and the pipe spacing (2 r^). The 
pipe length is determined from the fixed volume requirement: 

Lj - (3.58 X 10® ft®)/(r? - r®) mr (49) 
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The Biot number is given by — 

Ngi = hr^/kg (50) 

so that an expression for the steam-side heat-transfer coefficient, h, must be 
derived. In this example, we will first assume fully turbulent flow in the 
steam pipes that allows the use of the Nusselt equation — 

= 0.023 Npr®*^= (51) 

to derive the expression for h. Solving Equation 51 for h yields — 

h = 0.02 Npj.®*^ k (:^)®‘® n~®-® r^"^*® (52) 

where n = number of pipes with radius r^. Substituting in Equation 52 the 
appropriate properties for steam at 725“F (385®C) and 1000 psi, and setting 
rjj = 1.25 in. and m = 5 X 10® Ibm/hr yields — 

h = 5968/n®-® Btu/hr-ft^-®F (53) 

It can be seen from Equation 53 that and r^ are both dependent on the number 
of steam pipes to be used (Ngj^ “ n"®»® , r^ fixed). The pipe length, Lj, can 
be determined for any number of pipes by selecting n and calculating h from 
Equation 53, determining and using it to find rj from the Megerlln plots, 
then substituting rj and n into Equation 49 and solving for L^. 

Now that we have calculated a pipe length consistent with a given number 
of tubes, n, that are sufficient to remove the heat from the salt volume, we 
must consider if sufficient heat-transfer surface is present to heat the steam 
to the required outlet temperature. 

If it is assumed that the PCM remains at its melting point, the enthalpy 
change in the steam as it progresses through the pipe is described by — 

Rearranging Equation 54, integrating, and evaluating, given that at 
L = 0, t - 550*F (288®C) and at L - L 2 , t - 900“F (482"C) leads to - 

, _ (2.255) m Co 

^2 niHa (55) 
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The heat-transfer requirements will be satisfied when n is selected so that 
Lj L 2 > which meets both the fixed-volume and steam-temperature requirements. 

The solution will be obtained graphically here using Equations 49, 53, and 
55, and the Megerlin plots in Figure 45 to derive Table 15. 


B 


Nph-2. %o“15 

.9, tQ « 1.25 in. 

, AND m » 5 

X 10® 

Ib/hr, USING LIKCO3 AS THE PCM 


n 

h, Btu/hr-ft^-®F 

«Bi 

r Ll In 

• Ll » ft 

ta, ft 


3,000 

9.87 

0.79 

3.56 

49.3 

38.4 

1 

5,000 

6.56 

0.53 

3.25 

36.5 

34.7 

-3 

A 

J 

7,000 

5.01 

0.40 

3.00 

31.6 

32.4 


10,000 

3.77 

0.30 

2.75 

27.4 

30.2 


15,000 

2.72 

0.22 

2.50 

23.4 

27.9 

-i 

20,000 

2.16 

0.17 

2.31 

21.8 

26.3 


The values of Li 

and La are shown 

as a function of 

n in Figure 46. 

The inter- 

i 


section of the two curves, representing the design solution, occurs at n > 6000 
and Li > L 2 = 33.5 ft. Substituting n ■■ 6000 into Equations 49, 53, and 55 
yields h “ 5.67 Btu/hr-ft*-®F, ri ■ 3.13 in., and Li = La 33.4 ft. The design 
conditions for the TES operation can be met by solidifying 5 million lb of 
LIKCO 3 cyllndricaliy around 6000 2.5-in.-dlam by 33.4-ft-long pipes (spaced 
6.25 in. apart) over a 6-hr period. As a check, the total area open to steam 
flow is — 


A ■ nxr 

The linear steam velocity is — 

»5 


(6000) (X) (ij|^) 


204.5 ft 


(56) 


(5 X 10^ Ib/hr) (hr/3600 s) (0.6285 ftVlb) . ^ 

204.5 ft^ ^ 


(57) 


where the specific volume of steam at 100 psl, 725"F (385*C) is 0.6285 ft’/lb. 
The Reynolds number is — 


N' 


Re 


. dv£ . (2.5 in./12 in./ft)(0.43 ft/s)(0.6285 ftVlb)"V36C0 s/hr) _ 


0.0575 Ib/hr-ft 


which indicates fully turbulent flow and verifies the use of the Nusselt 
equation to derive the steam-side heat-transfer coefficient. 
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Figure 46. GRAPHIC SOLUTION FOR NUMBER OF PIPES AMD THEIR LENGTH 
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Although the above provides the design of. a TES system capable of meeting 
the performcnce requirements. It Is not necessarily the optimum system design. 
The use of different steam pipe sizes and the effects of conduction enhance- 
ment on the solidification process must be studied before an optimum performance 
design could be established. Such optimization must also include the determl- 

'r 

nation of the cost-effectiveness of the possible designs, but Is at present 
beyond the scope of this project. The example presented here does, however, 
show the Importance of the heat-transfer model In determining key parameters 
for any latent-heat TES system. 
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Task 5. Conclusions and Recommendations 
Conclusions 

1. Salts 

• Carbonates appear most attractive from thermal, chemical, and physical 
properties; LlKCOs selected as model system. 

• Chlorides appear most attractive from cost considerations. 

2. Heat-Transfer Modeling 

• Megerlln's solution gives good agreement between predicted 'and observed 
Interfaclal location during heat removal phase from PCM. 

• Model can be used for determining key design parameters for large-scale 
systems. 

3. Inert Additives to Control Volume Change (LlAlOg) 

• Small volume change observed with model system does not require their use. 

• Sedimentation does present a problem at concentrations of 20 volume percent, 
e Performance is affected; heat flux lower. 

• Inerts Increase system cost substantially while lowering specific capacity. 

4. Conduction Enhancement 

e Additions of stainless steel screen, fins, and wool, and aluminum wool 
(3-14 volume percent) Increase performance by 25%-32Z. 

e Wools appear most attractive of the enhancement concepts tested, particu- 
larly aluminum wool. 

5. Containment 

e Low corrosion of 304 and 316 SS with carbonates. 

• Carbonates show no creepage, little vapor loss. 

• Aluminum appears stable In the carbonate environment. 

6. Engineering Performance 

a Cycling easily accomplished within a 24-hour period, 
a Stable performance through 1400 hours of operation. 

OF POOR 

a Cycling has no apparent Influence on the salt properties. 
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• Performance unaffected by degree of superheating. 

• Performance was Improved 30% by the addition of conductivity enhancement 
in the form of eight fins. 

• Conduction enhancement provides more efficient utilization of heat- 
transfer area. 

• Established baseline average discharge heat fluxes for comparative testing: 

No Enhancement: 4000 Btu/hr-ft^ 

With Enhancement by Fins: 5300 Btu/hr-ft^ 

• Strong convective mixing effects Identified. 

• Low coefficient of expansion and related volume change on fusion of 
LIKCO 3 verified. 

• Strong dependence of heat flux on solid-salt thickness. 

• Heat fluxes as high as 40,000 Btu/hr-ft^ observed, representing an upper 
performance limit. 

7. Overall Conclusion 

• Alkali carbonates are suitable for hlgh-temperature TES applications 
from a performance, thermal, and chemical behavior standpoint. 

• Economic feasibility can now be studied. 

Recommendations for Further Research 

1. Salts 

• Study of lower-cost salts (commercial grade) 

• Wider range of melting temperatures 

e Incongruent compositions 

2. Modeling Studies 

c More detailed study of superheating and convective mixing effects on 
solidification and melting behavior. 

• Model verification with different sizes and shapes of containers 

e Model verification using a higher-temperature working fluid (higher Npi^) 

e Model verification using an alternative PCM 
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3. Inert Volume-Change Control Additives 

• No further study recommended at this time for carbonates 

4. Conductltm Enhancement 

e Further Investigation of SS or a3,umlinim in various forms, primarily 
wool or fixed-structure configurations 

• Corrosion testing of promising materials in different carbonate mixtures 

5. Containment 

• ' Long-term corrosion testing of current .containment materials in several 

salt mixtures under cycling conditions 

• Identification and/or development of lower-cost containment materials 

• Compatibility testing under thermal-cycling conditions of other candidate 
materials 

6. Performance Testing 

e Use of hlgher-tenq;>erature working fluids 

• Alternate passive heat exchanger configurations (multiple type) to 
establish performance in a more "commercial" system design 

• Larger-scale testing cf systems designed from the heat-transfer model 
e Long-term endurance testing (>10,000 hr) 

e Long- and short-term testing with alternate PCM'e 

t Long- and short-term testing with most promising conductivity promoters 

7. TES System Development 

s System studies to determine design constraint to be placed on systems in 
various uses 

s Cost and energy-saving analyses for promising concepts and economic 
feasibility 
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NOMENCLATORE 

heat-transfer area, ft^, in.^ 
constant (Equation 14) 
characteristic thicky.ess, ft, in. 

(half- :hickness of a slab, radiu-s of a cylinder or a sphere) 

heat capacity, Btu/lbm-“F 

3 3 

differential operator = v- + v , — 

oT o X 

wall thickness, ft 
tube diameter, ft, in- 

volume fraction of container occupied by molten salt 
mass velocity, ibiu'hr-ft^ 
acceleration due to gravity; g » 9.807 m/s^ 
surface heat- transfer coefficient, Btu/hr-ft^-®F 
heat of fusion, Btu/lbm 

V=i 

thermal conductivity, Btu/hr-ft-°F 
wave number, ft”" 
significant dimension, ft 

pipe length to meet fixed volume requirement, ft 
pipe length to meet energy transfer requirement 
dimensionless frequency; 

total mass, Ibm 
mass flow rate, Ibm/hr 



n = number of heat exchange pipes 

Np^ ® Fourier number (dimensionless time) ; %o ■ ar/a* 

Npp ■ dimensionless period; Np^ = aTg/a^ 

%u * Nusselt number (ratio of surface conductance to fluid conductance); 

%u • 
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Nph = phase-change number (ratio of heat of fusion to the sensible heat 
between the melting point and a reference temperature); 

^ph ” ^^f^i^^ps^^m “ ^s 

Np^ ■■ Prandtl number (ratio of kinematic viscosity to thermal diffusivlty); 
Npr “ v/a 

^Ra " Rayleigh number (ratio of buoyant to thermal forces) 

Np - Reynolds number (ratio of inertial to viscous forces); 

N,^ = vd/v 

Nbi “ Biot number (ratio of surface conductance to interior solid conduc- 
tance) ; Ngj ■ ha/kg or Ua/kg 

p ” perimeter, ft, in. 

R = r/r^j, radius ratio 

r •= radial coordinate, ft 

r* = radial Integration variable 

Q = quantity of heat, Btu 

q = heat flux, Btu/hr-ft^ 

t = temperature, “C, 

At = temperature gradient, ®C, *F 

T = absolute temperature, “R 

U = overall heat-> transfer coefficient, Btu/hr-f t*-'F 

V * fluid volume, ft"* 

V = velocity, ft/hr 

X ■ moving length coordinate, ft 

X ■ length coordinate, ft 

z • longitudinal coordinate, direction of fluid flow, ft 

Greek Letters 

a ■ thermal diffusivlty, ft*/hr; a ■ k/CpP 

8 ■ coefficient of volume expansion, *R“^; 8 "(Pj*-P)/pAT 

5 "a finite difference operator 
A ■ same as 6 
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X = wavelength, ft 

y = dynamic viscosity, Ibf-hr/ft^, cP 

V = kinematic viscosity, ft^/hr 

frequency, cps 

p = density, Ibm/ft* 

po = density of PCM, Ibm/ft^ 

p^ = density of liquid salt far away from heat-transfer surface 
T = time 

Tq = period 

(I) = circular frequency, rad/s; 0 )= 2irv 


Subscripts 


a 

= 

ambient 

c 

- 

container 

app 

s 

apparent 

f 

= 

fluid 

in 

- 

initial 

ins 

s 

insulation 

i 

= 

liquid (phase) 

m 

= 

melting 

s 

- 

solid (phase) 

surf 

= 

surface 

w 

s: 

wall 


Superscripts 

* per unit time 
— = average quantity 
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